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SECTIOH 1 
INTRODUCTION 


This report presents the results of an analysis effort performed to 
demonstrate the feasibility of employing approximate dynamical models and 
frequency-shaped cost functional control law design techniques for helicop- 
ter vibration suppression. Both fixed-gain and adaptive concrol designs 
based on linear second-order dynamical models were implemented in a detailed 
Rotor Systems Research Aircraft (RSRA) simulation to validate these active 
vibration suppression control laws. Approximate models of fuselage 
flexibility were included in the RSRA simulation in order to more accurately 
characterize the structural dynamics. The results for both the fixed-gain 
and adaptive approaches are promising and provide a foundation for pursuing 
further validation in more extensive simulation studies and in wind tunnel 
and/or flight tests. 


1 . 1 BACKGROUND 


One of the Important improvements which will influence a vilde accept- 
ance of rotorcraft for both commercial and military applications i3 the 
reduction of the mechanical vibration level. Reducing the vibration level 
will effect ride quality both directly and indirectly through its affect on 
noise level, a3 well as providing a more stable weapon delivery platform, 
and may reduce maintenance, cost, and weight. 

There are several significant sources of vibration in rotorcraft in- 
cluding rotors, transmission systems, and engines. The dominant source of 
vibration for most helicopters in forward flight, however, is the periodic 
loading of the rotor blades. When a helicopter transitions from hover to 
forward flight, it loses polar symmetry of the airflow through the rotor 
disc. The relative velocity of the air which the blades encounter is (to 


1 

i 


first order) a sinusoidal function of the blade azimuth with respect to the j 

direction of flight. These differences are further augmented by a helically < 

shed trailing vortex system, resulting in a rough ride for the blades and 
consequently the aircraft, even in smooth air. Vibration at high forward 
speed is further accentuated by compressibility and stall effects. 

Basically, the vibrating loads are periodic with the majority of the 
power concentrated around frequencies which are harmonics of n, the rotor 
angular velocity (h p ) times the number of blades (N) . 

0 - ( 1 . 1 ) 


Though vibration can occur at frequencies other than n, 20, 30, etc., due to 
structural resonances in the airframe and the blades, the disturbance power 
at frequencies other than harmonics of 0 is significantly below that near 0 
and its harmonics. This results in significantly lower vibration levels. 
Near 0 and its harmonics, however, rotor ana fuselage flexible modes can 
transmit the rotor hub load3 to various locations in the airframe with sig- 
nificant increases in amplitude and alterations in phase. 

Uncontrolled vibration levels in rotorcraft can exceed 0.25 g, a level 
significantly higher than that for fixed-wing aircraft. A example of uncon- 
trolled vibration for a MBB BO 105 helicopter in steady-state forward flight 
at 70 knots and during some maneuvers is show i in Figure 1-1 (taken from 
reference [1]). Notably, the accelerometer and rate gyro outputs have been 
pre-filtered with a three-pole analog filter with cutoff frequency 16 I'.z, 
yet the vibration level remains quite high. Current methods for reducing 
these vibration levels involve the use of mass-spring-dampei 1 systems for 
isolating and/or absorbing vibrational energy. These passive control 
systems are heavy, and are only effective at a few operating conditions. 
This has led to considerable research into 'active' techniques, where the 
rotor blades are excited at harmonics of 0 in order to reduce the 
aerodynamic loads. A summary of the findings of this research i3 given 
next. 


















1.2 ACTIVE CONTROL OF ROTORCRAFT VIBRATION 


Moat of the research to dato on active control of rotorcraft vibration 
has focused on multicyclic (or higher harmonic) control (N/rev and lt3 
harmonics) since the pioneering work of McCloud and Kretz [2], An excellent 
review of the state-of-the-art in multicyclic control as of 1982 wa3 given 
by Johnson C 33 - He characterized multicyclic control as that with: (1) a 
linear, qua3l-static, frequency domain model of the helicopter vibration 
response to control; (2) identification of the helicopter vibration model by 
least-squared error or Kalman-f liter methods; and (3) a minimum variance or 
quadratic performance function controller. Johnson also presented a 
tutorial on such control algorithms that use recursive parameter identifica- 
tion with linear feedback, algorithms which are generally termed self-tuning 
regulators. Discussions of other feedback strategies and helicopter vibra- 
tion in general can bo found in Johnson's text [H], 

The linear, quasi-static, frequency domain model of helicopter vibra- 
tion referred to requires further elaboration in order to accurately express 
the differences between this model and the linear 3econd-crder dynamical 
model employed in this analysis. Thi3 linear qua3l-static model, referred 
to in the literature as the 'T-matrix' model, is quite simply a linear 
static model cf the input-output relationship expressing the observed ac- 
celerations (z) as constant linear combinations of the Inputs (u): 

z - Tu . (1.2) 

This is clearly not a valid model over the entire frequency range of 
helicopter operation, but can be used ho accurately express the input-output 
characteristics at a 3lngle frequency. As discussed by Gupta in [5], 
restriction or an arbitrary transfer function description of the input- 
output relationship results in such a static relationship (enlarging the the 
field of scalars to include complex numbers). 

Since the static model is unknown a priori, it must be identified, and 
since helicopter dynamics are inherently nonlinear, the identification must 
be performed adaptively (on-line) or at a prespeoified number of flight con- 
ditions within the flight envelope. In order to restrict the identification 





to a specific frequency (or frequencies), the data input to the various 
identification algorithms arc effectively bandpass filtered by transforming 
into the frequency domain (FFT) and retaining only the responses (magnitudes 
and phases) at the frequencies of Interest. Thus, to state expllcity what 
is meant by the terminology "linear, quasi-static, frequency-domain model", 
the reference is to a static linear model relating the inputs and outputs 
which is designed to be valid only at a specific frequency (or frequencies). 

In this analysis, the linear static model has been extended to Include 
second-order dynamics between the inputs and outputs. This broader class of 
dynamical models is used in order to account for input lags due to blade 
dynamics as well as the effects of fuselage flexibility over a finite 
bandwidth near n (N/rev). Coupled with an appropriate control strategy, 
this class of models will admit identification and control of non-minimum 
phase systems (systems with transfer function zeroes in the right half- 
piano). Non-minimum phase systems are commonly found in the relationships 
between control surface deflections and achieved accelerations on 
rotorcraft, aircraft and missiles, 

Both linear time-invariant controllers and adaptive controllers were 
studied. The application of 'optimal* regulator theory to the helicoptor 
vibration control problem was successful primarily due to the utilization of 
frequency-3haped cost functions in the linear-quadratic regulator (LQR) 
design procedure [6.7]. These cost functions lead to control feedback whose 
energy is concentrated in the frequency band or band3 of interest. Time- 
invariant controllers were tested on the RSRA simulation both with and 
without fuselage flexibility effects included [8], As expected, fuselage 
flexibility radically altered the input-output dynamical relationships, in- 
dicating the need for considering these effects in any control design. 
Furthermore, the input-output characteristics (of the nonlinear RSRA 
simulation) were shown to vary significantly with flight condition. Though 
gain-scheduling techniques can be used to update the controller a3 the 
flight conditions change, other considerations such as configuration change 
which alter the input-output characteristics make investigation into the U3e 
of adaptive control laws potentially useful. 

The research reported here and by Gupta [5], differs from previous ap- 
proaches in two principal ways — in- the model form and in the control cost 



functional formulation and resulting mechanization. Previous work has con- 
centrated on static models in which the inputs and outputs are directly 
related by a coefficient matrix (T) whose parameters (elements of T) are 
modeled as stochastic random-walk processes. As outputs of an identifica- 
tion algorithm, these parameters may slowly vary with time. In this 
analysis, the input-output model is augmented to include second-order 
dynamics. Thi3 more complex model more accurately characterizes the 
helicopter's input-output relationships near M/rev. Using frequency 
weighted cost functionals in the estimation and control law synthesis places 
increased emphasis on errors in the frequency bands of interest. 


1.3 SUMMARY OF RESULTS 

Fixed-gain and adaptive control laws were designed for suppression of 
the narrowband vibration (at the pilot location) resulting frers periodic 
loading of the rotor blades in a simulation of the rotor systems research 
aircraft (RSRA). The particular version of the simulation used did not in- 
clude blade dynamics or aeroelastic effects, however the simulation wa3 
augmented with a model of eight flexible modes of the fuselage near n 
(H/rov). The results of the fixed-gain control design indicate a sig- 

nificant 'robustness' with respect to changes in flight conditions including 
speed, altitude, and weight (disc loading). _ 

The adaptive controller also performed remarkably well at all flight 
conditions, including a 35 second transition from 80 knots to 125 knots. In 
general, the convergence times for the adaptive controller were larger than 
those for the fixed-gain controller, but the steady-state RSS (root-sum- 
squares) acceleration at the pilot location wa3 insignificantly different in 
most cases. Further comparisons of results for the two strategies are dif- 
ficult to Justify in light of the fact that different control :onf iguration3 
were used for the two cases (four controls including the t.all collective 
were used in tne fixed-gain design while the tail collective was eliminated 
in the adaptive controller), and more fundamentally gince the solutions to 
the' vibration control problem used two different sets of underlying 
assumptions. However, it 13 interesting to note that in most cases the 





However, it is interesting to note that in most cases the adaptive control- 
ler required less control authority to accomplish the same reduction in 
vibration level at the 3ame flight conditions- 

The absence of a detailed rotor model including aeroelastlc effects may 
have conributed strongly to the similarity in performance of the two algo- 
rithms. It is likely that the inclusion of higher order aeroelastlc effects 
would have resulted in more pronounced model variations as a function of 
flight condition than were realized in the RSRA simulation. The more 
pronounced model variations would have led to Increased performance degrada- 
tion of the fixed-gain controller away from the design point. However, due 
to the periodic nature of the disturbances, the approximate model of the 
vibration dynamics is adequate for periodic disturbance suppression in spite 
of the underlying complexity (nonlinearity) of the actual helicopter vibra- 
tion dynamics. 


1.U SUMMARY OF REPORT 


The report is organized as follows. Section 2 describes the addition 
of fuselage flexibility equations to the GEHHEL. simulation of RSRA. Section 
3 discusses the fixed-gain control law design, it3 performance, and presents 
some preliminary results ■ on robustness of the control law with respect to 
variations in flight conditions away from the design point. A discussion of 
frequency shaped co3t functions is also included in Section 3. Section H 
describes the adaptive control algorithm in detail, and presents the results 
for various flight conditions. Section 5 summarizes the results and key 
features of the adaptive and fixed-gain algorithms with performance com- 
parisons where appropriate. Section 6 contains recommendations for future 
work. 
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SECTION 2 

MODELING RSRA ROTOR CRAFT VI B RATIO H DYNAMICS 


A previous study ([8]) and this current research have used a detailed 
Rotor Systems Research Aircraft (RSRA) simulation (GENHEL) to validate the 
effectiveness of frequency-shaped cost functional control laws for helicop- 
ter vibration suppression. Since flexibility effects can dramatically 
change the dynamical relationship between control inputs and observed vibra- 
tional accelerations, several fuselage flexible modes near n (N/rev) were 
added to the GENHEL simulation of the RSRA. While the GENHEL rotor model 
was not as complete as it might have been, the basic rigid body rotorcraft 
simulation with fuselage flexibility calculated accelerations similar in 
magnitude to those observed on the actual rotorcraft. Unmodeled actuator 
and blade dynamics coupled with nonlinear aeroelastic effects certainly 
would result in different magnitude and phase (most significantly transmis- 
sion zeroes) relationships between the inputs and the outputs of the 
simulation model from those of the actual helicopter in the narrow band of 
frequencies around n and its harmonics. However, the dynamical model used 
in the control designs has sufficient degrees of freedom (pole3 and zeroes) 
to accurately characterize the input-output relationships in the frequency 
intervals of interest for the purpose of periodic disturbance rejection! 


2.1 THE NASTRAN MODEL AND VIBRATION EQUATIONS 


The GENHEL helicopter simulation program assumes the fuselage is rigid, 
ignoring the effects of fuselage flexibility on vibration. In order to in- 
clude theso effects, a finite element modeling program (NASTRAN) wa3 used to 
estimate the natural frequencies and mode shapes in a finite element model 
of the helicopter. An appropriate subset of the modes was chosen for inclu- 
sion into the GENHEL program, effectively replacing the rigid body with an 
appropriate flexible fuselage. This section provides a mathematical 
description of the flexible fuselage and its inclusion into the GENHEL 
program. 


blank not filmed 
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Fop a set of dynamical equations of the form: 

MX + Kx - 0 (2.1 ) 

where M is a symmetric mas3 matrix, and K is a symmetric stiffness matrix, 
the NASTRAN program computes the natural frequencies (u's), a generalized 
stiffness matrix K, a diagonal generalized mass matrix [H] , and the 
generalized eigenvector matrices ( 4 >'s), such that: 

$ T K<> - K, (2.2) 

- H. (2.3) 

NASA Ames provided a NASTRAN model for the flexible fuselage of the RSRA 
helicopter. For the first 3H (lowest frequency) vibrational modes, natural 
frequencies, generalized stiffnesses, generalized masses and mode shapes 
(generalized eigenvectors) at various locations throughout the fuselage were 
calculated. The results are summarized in Appendix A. 

The effects of external forces applied at various locations on the 
fusolage can bo considered by modifying equation (2.1) to include a non-zero 
forcing term: 

MX + Kx - TF, (2.4) 

where T is a force distribution matrix, and F is a vector of external 
forces (in the x, y, and z directions in a body-fixed coordinate system) and 
moments (about the x, y, and z axes denoted L, M, and H respectively): 

F T - [ F , F , F , L, M, N ] . 
x’ y’ z’ ’ ’ 

Note that since the z-axis is coincident with the rotor shaft, no yaw mo- 
ments (neglecting bearing friction effects) are transmitted to the fuselage. 
Using equations (2.2) and (2.3), equation (2.4) can bo written in modal 
form: 

firi + n 2 n - <i> T TF f (2.5) 


10 - 



whero 


x - <j>n. (2.6) 

Since M la a positive, diagonal matrix, its Inverse is well-defined and the 
differential eouatlona can be written: 

h + M _1 n 2 n - m'Vtf, (2.7) 

where 

M _1 n 2 - diag{w 2 } (2.8) 

la a diagonal matrix of the natural frequencies. For the RSRA helicopter, 
theso modal differential equations can be written: 

T 
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where T la a rectangular (34 x 6) matrix of zeroes with an imbedded 6x6 

T 

identity matrix which picks out of <{> the vectors associated with the 

modal deflections at the hub. Thu3, defining the columns of T as T_, T„, 

x y 

Tp, T^, T m , , we can identify, for example, Tp $ as the NASTRAN state at 
z t x 

the hub, where the x-dlrecilon force is applied. Similar identifications 

hold for the other five degrees of freedom. Using the fact that the mass 

and stiffness matrices in equation 2.9 are diagonal, the equations decouple 

yielding: 



where 


"m * "m n m 


( 2 . 10 ) 


• ♦Ih F/f V 

♦oh ' [♦.h' F x>- WV- WV- W L1, W H >- 

and F la applied at the rotor hub. The vector $ is the vector of modal 
deflections at the hub for mode m. 

The outputs of interest are translational accelerations (specific 
forces) at the pilot's location. In terms of the vibration states at a par- 
ticular location on the fuselage, these outputs can be written as follows: 

y p - HX p . (2.11) 

where x p is the displacement of the fuselage at point p (eg. the pilot 
location) as a function of tine, and K p is the resulting acceleration. 
Those measurement equations can be written in terms of the modal states as 
follows: 


HC *1p *2p V * ’ V* ’ 


( 2 . 12 ) 


where y are acceleration measurements at the pilot location and a, is a 

P th 10 

column vector whose six elements are the deflections of the 1 mode at the 

pilot location for each of the six degrees of freedom. The measurement dis- 
tribution matrix H specifies the contribution of each degree of freedom to 
the accelerometer output. For an accelerometer mounted along the x- 
directlon, the H-matrix would be proportional to: 


H x *[ 1 0 0 0 0 0 ] 


Up to this point, the effects of structural damping have been ignored. 
In order to include these effects, a term proportional to n is added to the 
modal equations. The constant of proportionality is in turn proportional to 
the damping coefficient. Since typical structural damping ratios are be- 
tween 3 1 and 5$, the modal equations were augmented with a damping term 
giving 5? damping. These equations were discretized using zero-order hold 



equivalence and a sample Interval equal to the GENHEL simulation time 3tep. 
The resulting discrete equation is of the following form: 


'k+ 1 ,1 


a 1 i^k.I ' a 2»Vl,l * "olVl * b tl u k * b 2lVl ! 


where k is a time index, 1 Is a mode index, u k 13 the applied force at the 
hub at time k, and y . is the contribution to the accelerometer output at 

K p 1 

the pilot location at time k from mode i. The coefficients a 1 i» a 2i ,b oi ,b i i ’ 
and b 2i are the polynomial coefficients in the discrete-time z-transform 
transfer function representation of the zero-order hold equivalent input- 
output relation (rotor hub forces to pilot location accelerometer outputs). 
This equation was coded and added to the GEfJKEL simulation program. 


2.2 MODAL CONTRIBUTION AND MODE SELECTION 


The NASTRAN finite element modeling program was used to calculate the 
first 3H (lowest frequency) nodes for a finite element model of the fuselage 
of the RSRA helicopter. However, not all of these nodes contribute substan- 
tially to the vibration at the pilot's location. In order to significantly 
reduce the computational load, only a small subset of the modes which effec- 
tively contribute to the total vibration was chosen for inclusion into the 
GENHEL simulation. 

To select the modes with the greatest effect on the vibration at the 
pilot location, the amplitude of the response of each mode at the pilot 
location was estimated. The eight (8) modes with the largest responses near 
ft were Included. The modal responses were estimated by assuming that the 
vibrational forces in all six directions (three forces and three moments) 
were of equal magnitude and phase at ft where for the five-bladed RSRA;: 

ft - Nft^ - 110.7 rad/sec, 

N is the number of blades, five, and ft^ is the rotor angular velocity. 
Then, using reciprocity, the amount of excitation of each mode due to the 


13 “ 



six forces and moments was estimated by the 3um of the modal responses at 
the hub: 


6 

l 

, _ k^L 

mo “ r. r , 2 n 2.2 

M [(id - fl ) 
m m 


W k > 

ra 


(2.13) 


The 3um in the numerator is over the six degrees of freedom. The steady- 
state response in each direction (in the force and moment 6-space) was then 
easily calculated using the corresponding modal amplitude at the pilot loca- 
tion (d'jjjpC * > > * 


A 


mp 


■'tup^mr * 


(2.14) 


L U 

A (k), the k element of A , is a quantitative measure of the amplitude 
tap mp ^ 

of the vibration response at the pilot location due to equal forces in all 

directions at the hub. The sum over k of the absolute values of A (k) 

mp 

(i.e. the L,-norm of A ) was used as a cumulative measure of the vibration 

mp 

contribution from each mode. A summary of these values for the 34 NASTRAH' 
modes is presented in Tables 2-1 and 2-2. 

In order to minimize the computational load in the controller implemen- 
tations, only significant modes near £1 were included in the GENHEL. 
simulation. This resulted lr the majority of th-' vibrational energy in the 
outputs of the simulation appearing at n with greater than lOdB of attenua- 
tion at the higher harmonics. This, in turn, allowed the controllers to 
concentrate soley on vibration suppression at a single frequency G, thus 
minimizing the dimensionality of the problem to be solved and hence the com- 
putational load. The six modal displacements for the eight selected modes 
at both the rotor hub and pilot location are given in Table 2-3. 
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TABLE 2-1. 


ESTIMATED MODAL RESPONSE A_ AT THE PILOT LOCATION 

P 


Mode 

(J 

No. 

(rad/s) 


Relative Modal Deflection at Pilot Location 


0.000 

0. 17531-05 

-0.11488-14 

0.1**38-12 

0.1 C 5°E-17 

C.::4J8-15 

0. 

0.005 

.0. 36138-14 

0.17538-03 

0.146BE-13 

-0.44538-15 

0.47525-17 

0 . 

0.000 

-0. 77548-14 

-0.44198-15 

0. 19048-05 

0.17948-16 

0. <3828-15 

0. 

0.030 

-0.25698-04 

-0.43528-06 

0.19038-05 

-0.6404E-07 

-0.64738-10 

0. 

0.000 

-0.2340E-04 

-0.1E288-03 

0.1 1678-05 

-0.24138-09 

0.65938-03 

0. 

0.000 

0.19648-06 

-0.14898-05 

0.29738-06 

-0.9S97E-03 

0.71518-11 

0. 

27.407 

0.25938-06 

-0.13398-07 

-0.13998-05 

-0.1702E-08 

-0.1492E-07 

0. 

21.032 

-0.14728-06 

0.65 648-04 

-0.64198-07 

0.26618-03 

0.4340E-0? 

0. 

<8.0*4 

0.86498-06 

-0.93678-03 

-0.22368-05 

-0.7277E-08 

-0. 38338-07 

0. 

59.641 

-0.57138-07 

0.74548-04 

0.2132E-C7 

-0.1250E-V8 

0.5*258-10 

0. 

64.732 

*. 12028-06 

-0.29738-06 

-0.15718-06 

0.42748-03 

-0.12068-03 

0. 

Cl .099 

0.20148-03 

0.44268-06 

0.52338-05 

0.32078-07 

C.1B20E-06 

0. 

93.432 

0.17748-03 

-0.2333L-0S 

-0.17348-05 

0.75E4E-07 

-0.44338-09 

0. 

97.437 

-0.13018-03 

-C. 14368-05 

-0.1672E-C5 

C. 44318-07 

-C. 24918-03 

0. 

U0.?<o 

-0.45108-06 

0.79168-06 

0.117S8-0S 

-V.3932E-07 

-2.21008-07 

0. 

117.015 

-v. 36168-04 

0.59258-06 

0.1908E-05 

-0.60448-07 

0.42438-07 

0. 

120.364 

-0.16058-05 

C.3SSDZ-05 

0.9501E-05 

-0.47278-07 

0.11348-06 

0. 

126.153 

0-12428-05 

-0.15308-05 

-0.29368-05 

0.74378-07 

0.53208-07 

0. 

121. .'14 

-1.42418-06 

0.22078-05 

-0.95758-05 

0.21228-06 

-0.43448-07 

0 . 

140.716 

-0.25298-06 

-C. 20408-06 

-0.13968-05 

0.4316E-03 

0.73128-07 

0 . 

156. '12 

-0.74998-07 

-0.29498-08 

0.19008-04 

0.12558-03 

-0. 39978-09 

0. 

165.927 

0.56578-08 

0.10518-07 

0.36*78-07 

-0.36628-09 

0.5*928-10 

0 . 

147.302 

-0.1Gc't-07 

-0.16428-06 

-0. 29238-06 

0.15938-07 

-5.11868-07 

0 . 

169.942 

0.27468-06 

-0.24358-07 

-0.46968-06 

-0.11088-67 

0.7121E-08 

0 . 

132.112 

0.21578-03 

-0.26328-07 

0.36168-05 

0.12128-07 

0.33138-03 

0 . 

107.013 

0.20738-06 

0,18698-07 

0.89978-05 

-0.51788-09 

0.19678-07 

0 . 

191.659 

-0.19648-06 

0.53998-07 

0.12208-06 

-0.22498-07 

0.56808-09 

0 . 

205.185 

C. 21098-06 

0.36748-0', 

v.o0365-Oc 

-0. 10668-07 

0. 40495-0S 

0 . 

217.013 

-0.2072F- ' S 

' -0.23608-07 

0.36*48-0" 

0.32128-0B 

-0.22825-07 

0 . 

224.810 

0.11778-05 

-0.11468-05 

0.46778-05 

-0.13708-06 

0.15965-06 

0 . 

2cj.73? 

-0.14788-05 

0.39238-06 

0.3691E-05 

-0.11058-06 

-0.11368-07 

0 . 

237.447 

0.15058-06 

C.9353E-0? 

-0.12928-05 

0.57558-07 

-C. 73048-07 

0 . 

233.716 

0.3S36E-06 

C. 91778-07 

0 7429E-06 

-0.24318-07 

0.12338-07 

0 . 

243.707 

-0.137CE-C4 

-0.42468-04 

-0 ' ,688-06 

0.6102E-03 

C. 11935-07 

0 . 









TABLE 2-2. CUMULATIVE AMPLITUDE RESPONSE AT THE PILOT STATION R 

mp 


• Mode 
Ho. 

'J 

(rad/s) 


R 

rap 

1 

0.000 

0.17531-05 

0.17532-05 

2 

0.000 

0.17531-05 

0.17532-05 

3 

0.000 

0.18043-05 

0.1 8042-05 

4 

0.000 

0.2709: -05 

9.4662E-05 

5 

0.000 

0.1410: -05 

0.51212-06 

6 

0.000 

0.1 9631-05 

0.65192-06 

7 

27.607 

0.16”:-05 

0.94952-55 

8 

31.032 

0.873 18- 56 

0.56232-05 

? 

45.046 

0.31523-05 

0.11772-04 

10 

59.661 

0.37532-06 

0.32032-05 

11 

64.732 

0.67042-06 

0.1C842-G5 

12 

01.09? 

0.79092-05 

0.37171-04 

13 

93.452 

0. 44152-05 

0.22652-04 

14 

97.437 

0.48552-05 

0.32972-04 

15 

110.946 

0.24 1 l£-05 

0.21342-04 

16 

117.015 

0. 29662-05 

C.2101E-04 

17 

120.3:4 

C.ef 241-Cr 

C. 29422-04 

IB 

126.15E 

0.58422-05 

0.16572-04 

1? 

131.614 

0.12562-04 

C. 29872-04 

20 

140.716 

0.19232-05 

0.13472-04 

21 

156.612 

0.23322-06 

0.20162-05 

22 

165.927 

0.57-352-07 

0.11222-05 

23 

167.502 

0.59122-06 

0.56712-05 

24 

169.942 

0.7S57E-06 

0.26328-05 

25 

182.112 

0.40552-06 

0.35332-05 

Zi 

icr rn 

« *» - • >*H 

0. 11512-73. 

C.5342E-:: 

27 

191.659 

0.39S4E-06 

0.30402-05 

23 

209.185 

0.932?2-06 

0.83022-05 

2? 

217.018 

0.30462-06 

0.29612-05 

30 

224.810 

0.72962-05 

0.14902-04 

31 

225.737 

0.56832-05 

O.E374E-05 

32 

237.447 

0.16712-05 

0.70192-05 

33 

238.716 

0.12302-05 

0. 37282-05 

34 

243.707 

0.7472E-06 

0.1817E-05 


16 - 
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TABLE 2-3. MODAL DISPLACEMENTS OF SELECTED MODES AT HUB AND PILOT LOCATION 


Mode 

Displacement at Hub 

No. 

• F 

X 

F 

y 

F 

2 

L 

M 

N 

12 

0.344E400 

0.160E-C1 

0.254E-0! 

0.38EE-03 

-0.880E-C2 

G.333E-03 

13 

0 . 106E40C" 

-0.427E-01 

-C.126E-02 

-G.237E-02 

-0.314E-C2 

-0.496E-03 

14 

-0.231E400 

-0.389E-G2 

0.396E-01 

-0.374E-03 

0.617E-02 

-0.214E-03 

13 

0.490E-C1 

0.313E-C1 

o.ioc>e-g: 

0.2G3E-C2 

— 0.113E-02 

0.655E-04 

16 

-0.3S1E-01 

-0.673E-C1 

-0.62CE-02 

-0.2411-02 

0.140E-C2 

-0.V76E-03 

17 

-0.180E400 

0.S22E-01 

-0.571E-C1 

0.2B2E-02 

0.546E-C2 

-0.559E-02 

18 

-0.137E400 

0.28iE-01 

-0.990E-01 

0.C81E-03 

0.434E-02 

0.229E-02 

1? 

0.17EE-02 

0.464E+00 

0.124E-C2 

0.199E-01 

-0.39E.-04 

-0.104E-01 


Mode 

No. 


Displacement at 

Pilot Location 


F 

X 

F 

y 

F 

z 

L M 

N 

12 

-0.542E-01 

0.119E-01 

-0.1411. 

-0.843E-03 0.49CE-02 

O.000E400 

13 

-0.7S3E-01 

-O.123EIO0 

0.766E-01 

-0.335E-02 -0.196E-04 

C.C00E400 

14 

0.455E-01 

-0.435E-Ci 

0.56SE-01 

-0.134E-02 -0.755E-04 

O.OOOEK-O 

13 

0.211E-01 

0.371E-01 

-0.352E-01 

C.1E4E-02 -0.516E-C3 

O.OOOE4CO 

16 

0.172E-01 

0.283E-01 

-0.90UE-01 

0.2BCE-02 0.202E-0" 

O.oOOErCO 

17 

0.S16E-01 

0.132E400 

-0.323E-01 

0. J61E-02 0.383E-C-2 

O.OOOE-tCO 

18 

-0.730E-01 

-0.923E-01 

0.177E+00 

-0.461E-02 0.3ME-02 

O.OO0E4OC- 

19 

0.142E-01 

0.737E-C1 

0.321E400 

-0.1 03E -Cl -0.14K-C-2 

O.OOOl+OO 


17 - 





Figures 2-1 and 2-2 compare the rigid body vibration at the pilot loca- 
tion (a root-auo-square (RSS) of the three force components) with the same 
vector sun of accelerations at the pilot location with fuselage flexibility 
modeled with the eight 3elected modes. The erratic behavior during the 
first second is due to the transients in the GENHEL simulation. The states 
in the GENHEL model are not initially in steady-state. A comparison of 
these two figures clearly indicates the importance of fuselage flexibility 
in terms of its contribution to the vibration at the pilot location (at 
least for this simulation!). 



Figure 2-1. RSS Acceleration at the Pilot Location for Rigid Body Model 
(120kt3/100 ft) 



Figure 2-2. RSS Acceleration at the Pilot Location for Rigid Body with 
Eight Mode Fuselage Flexibility Model (120 kts/100 ft) 





■NT 



Figure 2-3 shows the output of the GENHEL simulation with the eight 
flexible modes included at a flight condition of 120 knot3 forward speed and 
an altltudo of 1000 feet. FFT's of the last three (3) seconds of the ac- 
celerometer outputs clearly manifest the concentration of vibrational energy 
at £1. The power at harmonics of 0 is more than 10dB below that at 0, as it 
was designed to be. Inclusion of fuselage flexible modes at or near har- 
monics of Cl would certainly have contributed to Increased power at those 
frequencies, resulting in a increase in the computational load in the con- 
trol implementations. 
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Figure 2-3= Acceleration Outputs and Power Spectra at Pilot Location with 
Eight Mode Fuselage Flexibility Model (120 kts/1000 ft) 
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SECTION 3 

FIXED-GAIH' CONTROLLER FOR RSRA ROTORCRAFT VIBRATION SUPPRESSION 


The basic steps involved In the design of the fixed-gain controller 
include linearization of the nonlinear model of the helicopter vibration 
dynamics about the desired operating point and measurement model lineariza- 
tion at the same operating point. Appropriate performance specifications 
are selected, a controller structure chosen, and parameters in that struc- 
ture such as estimator and controller gains calculated based on an 
optimization procedure. In modern control theory, the LQR (linear- 
quadratic^rogulator) design procedure leading to steady-state solutions of 
Riccatl type equations is a common choice for controller design. The mini- 
mization of a weighted squared error results in a linear state feedback 
control law with constant gains in steady-state. In estimation theory, the 
dual concept is linear-quadratlc-Gausslan (LQG) estimator design. Thus, the 
two major steps in fixed-gain controller design are: 

1) Generation of linearized dynamical and measurement models, and 

2) LQ design for estimator and controller feedback gains. 

There are several methods for generating linear models from nonlinear 
dynamical systems. The basic differences involve the choice of the ap- 
proximating model form . Utilizing most of the physics underlying the 
problem, complex 'model forms can be arrived at and the parameters in these 
models (stability derivatives) Identified from input-output measurements at 
desired operating points. This method was used in preliminary analysis and 
will be discussed in a little more detail in Section 3.1. However, such 
models usually Involve a large number of states and questions of observ- 
ability, controllability, and computability become salient Issues. A second 
approach involves identifying parameters in a low order input-output equiv- 
alent model form. This naturally reduces the required computational load 
since fewer parameters are involved, at the price of a more approximate 
dynamical model whose ability to predict future system outputs given the 
pa3t states and Inputs will certainly be inferior to that of the more com- 
plex higher order model forms. However, as alluded to earlier on several 
occasions, the helicopter vibration problem is one of controlling outputs at 



a few known frequencies. The approximate model of the helicopter vibration 
dynamics in therefore required to be accurate only at those frequencies. A 
globally valid model is not necessary as long as the control authority is 
limited to those frequencies a3 well. 


3.1 IDENTIFICATION OF LINEARIZED MODELS OF THE RSRA VIBRATION DYNAMICS 


This subsection presents an overview of the detailed linearization pro- 
cedure used as a first attempt at obtaining a linearized model of the RSRA 
vibration dynamics. With the addition of eight fuselage flexible modes, 
however, the detailed linearized modol became excessively large and sub- 
sequently a low order input-output equivalent model was identified. The 
details of thi3 procedure are also dlscu33ed. With, the modol parameters 
Identified, LGR control design was performed by choosing cost3 on thG 
frequency-shaping filter states and solving for the optimal gains as dis- 
cussed in the next subsection. The fixed-gain controller design procedure 
is outlined in Figure 3 - 1 . 

3.1.1 Stability Derivative Estimation U3lng the RSRA Simulation 

Based on the approach described in [2], linear models of the RSRA 
helicopter were constructed by computing stability derivatives of the 
helicopter rigid fuselage and rotor model at various operating points. 
Small perturbations in the states of the model which included states as- 
sociated with the six rigid body degrees of freedom and a flapping mode for 
each of the five blades resulted in variations in the forces and moments at 
the hub, leading directly to an estimate of the associated stability 
derivatives. The performance of controllers designed with these models was 
entirely satisfactory when applied to the nonlinear, r igi d body simulation 
of the helicopter. However, severe degradation in performance resulted when 
flexible modes were added to the simulation. In fact, the fixed-gain rigid 
body model controller design was destabilizing at several flight conditions 
when applied to the flexible helicopter model. 
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The destabilizing effect of Including the fuselage flexible modes was 
not unexpected, and demonstrated clearly the need for inclusion of these 
modes in the control design. To Include the effects of the flexible modes, 
each mode with a significant contribution to the vibration was added to the 
dynamical model, and the measurement models altered to include these 
flexible mode contributions. Since the model order increases by two with 
each additional mode, the inclusion of a significant number of modes 
resulted in the stability derivative calculation becoming computationally 
burdensome and the controller implementation virtually infeasible (the num- 
ber of computations at each time step for fixed-gain controllers being 
proportional to the square of the number of states). 


3.1.2 Identification of a Low Order Linearized Model of the RSRA Simulation 

To reduce the number of states in the linearized model of the RSRA, 
several approaches can be taken. The high order model including all the 
flexible modes of interest can be reduced by one of several procedures for 
eliminating combinations of states which are weakly observable with the 
measurements available and/or weakly controllable with the Inputs available. 
This presumes, of course, that the effort to acquire such a detailed model 
ha3 already been expended. However, based on the observation of the rotor 
frequency is fixed in most helicopters, an alternate procedure was chosen. 
Since most of the power in the vibrational forces is at this frequency and 
its harmonics, approximate second-order models adequately describe the 
input-output relationships required for controller design in the frequency 
band of interest. As discussed in Section 1, this class of models Is an ex- 
tension of the class of linear static models to include control lag effects 
a3 well as approximate impulse response, characteristics. - The parameters in 
these models must be identified from input-output data before a fixed-gain 
controller can be designed. 

In almost all helicopters the rotor angular velocity (CJ^) Is kept 
nearly constant. Thus, the vibration is concentrated in narrow frequency 
ranges, and second-order systems can be employed to adequately represent the 
transfer function gain and phase over each interval. As previously dis- 
cussed, the fuselage flexible modes are the most significant contributors to 
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tho vibration at the pilot location. Since the objective of the analysis 
was to demonstrate the feasibility of various conti ol designs, the problem 
was simplified to Include only major contributors around ft, resulting in in- 
significant power at higher harmonics. Thus, a single second-order model of 
the RSRA vibration dynamics, designed to be accurate over a finite frequency 
interval around ft, was chosen for each accelerometer measurement. 

To excite all the flexible modes around ft for the purposes of model 
parameter identification, control deflections containing frequencies in a 
band above and below ft (approximately^ 20S, i.e., 90 rad/sec < u < 130 
rad/sec) were input to the RSRA simulation. Tho resulting accelerations in 
the x, y, and z directions along with the four control inputs and the 
'known' disturbances, costj; and 3in^, wore then used as outputs and inputs 
respectively of a linear autoregressive model with exogenous lnput3 (ARX). 
Sine and cosine of the rotor azimuth angle were included in order to account 
for the periodic disturbance of unknown amplltuao and phase. Algorithms for 
identification of the parameters in ARX models assume that the 3tocha3tic 
processes entering Into the models is white and Gaussian (WGIJ). Colored 
(time-correlated) noloe such as the periodic disturbances which are the 
sources of tho vibration must be taken into account ( pre-whi toned) , other- 
wise biased parameter estimates result. In this formulation then, the 
variables In the ARX model of the input— output dynamics were; 


[a , a , a ] , 
x* y' z ’ 


CB 13* V A ls’ 9 tr 3, 


and 


T 

w 


[cos<!>, sin<»] . 


For these three measurements and four control inputs, a sixth order 
linear discrete ARX model with fully coupled modes may be written as: 


Vi ■ A 1 y k + A 2 Vi + B 1 u k * 3 2 Vi + r w k U-V 

where the A^'s, B^'s and T are matrices containing the parameters to be 
identified. The A^'s are 3x3 matrices, the B^'s are 3x4 matrices, and r 
i3 a 3 x 2 matrix. For future reference, by defining; 
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(3.2) 


T k " A 2 Vi + B 2 Vi 


equation 3.1 can be written in state-space form as follows: 


'k +1 


k +1 


“ A 1 

1 


" y k" 


r* 

B ! 


0 


T, 


B„ 

L 2 

_ 


k_ 


L 2 J 


u k f 


( 3 . 3 ) 


For the purpose of least-squares identification, equation 3.1 was expanded 
for the measurement a^ as follows: 

(a xVl " a 11 (a x>k + a 12 (a y>k + a 1 3 (a z ) k 

* a 11 (a xVl * a 12 ( Vk-1 + a 13 ^ a z \-1 

+b 11 (B ls ) k +b 12< 8 c )+b 13 (A 1 3 ) k +b l4 ( 9 tr ) k 

+ b li (B ls*k-1 + b 12 (0 cVl + b i3 (A l3^k-1 + b 1 4 (G tr 5 k-1 

+ r i1 cosij) k + r i2 3ini!) k . ( 3 - ^ ) 

where the a i j » a ' jj j .» b ]t j • and r ij are th3 elements of , A 5 , B 1 , and r 

respectively. A linear least squares solution of equation 3 . H was used to 
find the parameters from a data record of 2500 points ( 2.5 seconds of data 
sampled at 1000 points/second). The same technique was used for obtaining 
ARX model coefficients for the remaining outputs a and a . The dynamic 

y z 

matrices for the identified linear model for a steady-state level flight 
condition of 120 knot 3 forward speed at an altitude of 100 feet are given in 
Appendix B. 


3.2 FIXED-GAIN CONTROL LAW DESIGN 


The performance objective of the fixed-gain controller was to eliminate 
the vibration at the pilot location which for the purposes of this study was 
concentrated at fl (M/rev). Since most of the vibrational energy was con- 
centrated at a 3ingle frequency, the control design was based on a low order 
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linear model which designed to be accurate near the vibration frequency. 

The fixed-gain frequency-shaped controller employed utilizes a state-space 

modol of the plant and optimizes a cost function which places a large 

penalty on the outputs (accelerations at pilot location a , a , a ) at the 

x y z 

vibration frequency. 


3.2.1 Derivation of the 


Controller 


In order to place an increased cost on the outputs of a dynamical sys- 
tem over a specified frequency range, a modification to the standard LQR 
cost functional is required. In this subsection, the details of modifica- 
tions required are discussed along with the implementation in the overall 
controller design. The continuous-time analog of the ARX model of the 
helicopter discussed in the previous subsection can be written in state- 
space form as follows: 


k - Fx + Gu + Tw 


y » Hx + Du 


(3.6) 


where y is the measurement vector of accelerometer outputs, u is the 
vector of control inputs, and w is the narrowband disturbance vector (cf. 
equation 3*3). Subject to these dynamical constraints, minimization of the 
following frequency-shaped cost functional is desired. 


sr l.?' 


(Jw)A( jio)y(Jw) + u (Ju)B(Jw)u( ju)]du 


(3.7) 


where 


A(Jm) - P m (J“)P m (ju) , 


P m (Jw) - P(Ju)diag {a 1 , a 2 , ...}, 


h Q s 2 + h^c^s + h 2 n 2 
3 2 * 2c 2 ns + n 2 


( 3 . 10 ) 


// • 



lOC-MACMTTUCC Ri 


and 


B - I 


The'Oj's are scalar weighting factors and n is the vibration frequency. 

The location of the poles and the- zeros of the scalar frequency-shaped 
penalty function, P(s), are determined by h Q , h 1 , h 2 , Cj and ? 2 . The log- 
magnitude and phase plot of P(s) for two sets of h' s and c's are shown in 
Figures 3-2 and 3-3. As can be seen in Figure 3-2, both the magnitude and 
slopes of the phase plot of the penalty function at u » (i are Infinite 
(since ; 2 -o). These properties manifest a lack of robustness to modeling 
errors and a significant performance sensitivity to the choice of center 
frequency. Using a small amount of damping makes the system more robust to 
estimation and modeling errors at the expense of vibration reduction at the 
center frequency. As can bo aeon in Figure 3-3, the 3-dB bandwidth of the 
frequency response of the penalty function P (for c -0.5?) is 2c8. Note 
that setting h Q -h 2 «0 and hj-1 results in unity gain at u - 0. 


4 



Figure 3~2. 


Frequency Response of the Frequency Shaped Filter 
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Figure 3^3. 


Frequency Response of the Frequency Shaped Filter 

for ,2C£3 (h -h -0, h -1, c - C -0.5J) 

3 + 2c«s + ft 


In order to obtain a state-space model for the frequency shaping fil- 
ters, these equations are transformed to time domain. This is accomplished 
by defining: 
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/ /. 

/■; 


-r 


z, “Py^ » i M i > 2 , • « * i N 


(3.1D 


where M is tha number of outputs to be frequency shaped. Using (3.10) in 

(3.11), 


hgS^ + h^cfls + h 2 n £ 
5 . - - 
s ♦ 2?ns + a 


(3.12) 


In the time doraa? n, (3.12) becomes 


\ + 2C0 i i * n 2 z x - h 0 y t * 2;nn 1 * 1 + n 2 *!^ 


Equivalently, (3.13) can be written in state-space form: 


where 


*f 3i ■ F fs x rs i + Wi • 


Z i ’ H f3 X f3 1 * Vi ’ 


(3.13) 


(3. I'D 
(3.15) 


0 

“ 

i 


3 1 

„2 

-2cfl 

* G fs“ 

-a 

_ B 2_ 


F fs- 


B, - 2ca(h 1 - h Q ) 


Hfs“ Cl 03 


and 


»2 " h 0 * “* 2 < h 0 - h ,»' 


( 3 . 16 ) 


Finally, substitution of (3.9) and (3.11) in (3.7), and using Parseval*3 
Theorem, equation (3.6) can be rewritten: 


fta 

IT 2 2 T 

- j [z dlag (c^, a 2 , ...}z + u Bu] at . 


(3-17) 


: \ 


' ■ j 

>. 


ThU3, the frequency-weighted optimization problem is transformed into a 
standard linear-quadratic-regulator (LQR) design problem. 


pr 
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As discussed in previous sections, the control law design employed a 
large cost on the outputs at the vibration frequency. As discussed above, 
the frequency-shaped co3t functional can be rewritten in the standard LQR 
form by augmenting the state-space dynamical model with additional states 
corresponding to those of a state-space time-domain realization of the fre-, 
quency weighting function. The LQR design for the augmented system then 
considers cost on the frequency-shaping filter states and controls only. 

The frequency-shaping filter used in the fixed-gain control design was: 


H(s) 


s 2 ♦ 


(3.17) 


and is shown in Figure 3-1. Transformation into the discrete-time domain i3 
required and was accomplished using the bilinear transformation (trapezoidal 
integration): 


_ 2(z-1) 

3 “ T ( z+ 1 ) * 


(3.18) 


where T is the sampling interval. Substituting equation 3.18 into 3.17, 
the frequency-shaping filter in the discrete- domain can be written: 


H(z) - 


V * V * b 2 

2 

z + ajZ + a 2 


where: 


b Q - b 2 - [1 + ^ 2 ] , 


(3.19) 


and 


2b, 


T 2 fl 2 - 4 

? 2 ' 

t n + 4 


Transforming this transfer function into a state-3pace description is 
readily accomplished: 


7 


V •« 4 

f 


\ 


\ 

\ 

l 


-h 

V? 
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x(k + 1 ) - F x(k) + 3 u(k) + r w(k) 
. P r P 

yCk) - H p x(k) ♦ D p u(k) , 


.'3.21) 


where y 13 the measurement vector of accelerations at the pilot 
location. The subscript ' p' is used In equation 3.21 to denote 1 ilicopter 
vibration dynamics giant matrices which ... e hereafter assumed to be In 
the discrete domain. Appending frequency shaping filter states to the pi nt 
states yields the following combined state-3pace description of the system 
which can be used to design the control law. For the purposes of illustrat- 
ing the form of the system matrices, two raea3u-enent3 and two inputs are 
assumed. 





G. H . 0 

fs pi 


I 



I ° 

G fs H P 2 j 0 F fs 


G f3 D P 1 

G fs D P 2 



D 

, D - 

Pi 

° P 2 


and 

Y - 

y 



+ 

~D u " 
P 



z 


_ H fs x fs i _ 


- D rs y _ 


( 3 . 22 ) 


The measurement equation is not In standard from, but can be reorganized to 
yield: 


? 


H 0 


X 


n «n 

D 

P 


P 

+ 

P 

. D fs H p H fs. 


X 

I 


D. D 

L f3 pJ 


u . 


(3.23) 


The cost function now can be written in terms of Y's and u's, or X's and 
U'3, where X f - [Xp, x^] . 


J - l [Y r (k)A Y(k) * u T (k)Bu(k)] 
k-0 


(3-2U) 


where 


a y - 


o ' 0 

! a i 

0 | 1 a. 


(3.25) 


Equivalently, 

0» 

J - I [X T (k)A y X(k) ♦ X T (k)A Xu u(k) ♦ u T (k)Aj u X(k^ ♦ u T B„u] 


k-1 


u 

(3.26) 


where 


A y - K A y H , 


and 


A Xu" H V 


Because of the feedforward terms present in tho outputs Y, there is cross- 
coupling between state costs and control co3ts. Thl 3 is evident in the 







latter description. Instead of optimizing ( 3 . 26 ) , an equivalent cost func- 
tion (J) can be minimized which contains no cross-coupling between states 
and controls. 


where 


J 


l CX T (k)A Y X(k) ♦ u T (k)Bu(k)] 
k-1 X 


B - B ♦ D T A x D 
A x - A x " A^D B 'd Ajj 
F - F - G B _1 D A X H . 
The control gain is also modified to give: 
C » C ♦ b _1 d t a x h . 


( 3 . 27 ) 


A derivation of these modifications is presented In Appendix C. 

Using tho plant identified by the least-squares identification method 
described in tho previous section, optimal closed-loop regulator gains were 
obtained via solutions to the 3teady-state algebraic Riccatl equations 
resulting from the minimization of the cost function J. Cost matrices A^ 
and B u were determined 3uch that the closed-loop eigenvalues of the model 
had damping between 1? and 5$. This restriction on damping was empirically 
found to balance the controller performance in all three axes. Open-loop 
eigenvalues of the augmented system, cost matrices used, and closed-loop 
eigenvalues and feedback gains calculated for a 120 knot speed and 100 ft 
altitude flight condition are giver, in Appendix B. The performance of the 
controller at this design point 33 well as the performance at various flight 
conditions and disc loadings away from the design point are discussed in tho 
next subsection. 
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3.3 FIXED-GAIH CONTROLLER PERFORMANCE 

This section presents tho results of the fixed-gain controller design 
and implementation. The robustness of the design to modeling errors is also 
discussed. In order to maintain continuity in the discussions, the figures 
are presented at the end of the section. 

3.3.1 Controller Performance at the Design Point 

Having identified a simplified linear model of the RSRA vibration 
dynamics and designed a fixed-gain 'optimal' controller as discussed in the 
previous subsections, the GENHEL simulation was augmented with the feedback 
control design and the controller performance investigated. Figure 3-4 
presents time histories of the three accelerometer outputs and the four 
feedback command control inputs for four (4) seconds of simulated flight at 
120 knots and 100 feet altitude. Ho process or measurement noise was in- 
cluded in the fixed-gain controller simulation. 

During the first second of tho simulation, no control authority wa3 ox- 
erclsed in order to allow the simulation time to reach a steady-state 
vibration level in all three axes. The transients are clearly manifest in 
the accelerometer outputs in Figure 3-4. These accelerometer outputs were 
also filtered to remove the center-of-mass specific force components which 
were of little interest in the vibration control problem. Similarly, the 
control deflections shown do not include the control deflections required to 
trim the RSRA at the indicated flight condition. 

As expected, the controller performance is excellent at the design 
point as indicated by the rapid convergence and magnitude of vibration 
suppression. The vibration levels in all three axes are reduced to ap- 
proximately 10' of their original levels. The root-3um-3quare (RSS) 
acceleration output in Figure 3“4 indicates the overall performance of the 
controller. 

3.3.2 Control Law Robustness 

To investigate the robustness of the controller designed at the 120 
knot speed and 100 feet altitude flight condition, the same control gains 
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were used and the simulation run at different velocities, altitudes, and 
disc loading. Tho results are shown In Figures 3-5 through 3 — 1 4 . 

Figures 3-5 through 3-8 show the vibration reduction achieved at dif- 
ferent helicopter altitudes at the same 120 knot speed. Little degradation 
in performance of the 'controller 13 seen for altitudes below 5,000 ft. A 
significant degradation 13 noted at 7,000 feet altitude (Figure 3 - 8 ), but 
still the overall vibration level is reduced to about 25? of uncontrolled 
level. It should be noted that tho low frequency low amplitude modulation 
of tho measurement and control time histories shown in Figure 3-8 is not an 
instability in the controller, but rather an artifact due to the subsampling 
of the program outputs required to reduce the number of points to be 
plotted, and represents an excellent example of the phenomenon of ’aliasing’ 
in sampled data systems. 

As evidenced by the control time histories in Figures 3-4 through 3~8, 
the major effect of operating the controller at increasing distances from 
its design point for this RSRA simulation is a significant increase in the 
control authority required to achieve the vibration reduction. A probable 
cause for the increased control deflection requirements is the decrease in 
control effectiveness as the air density decreases. The air density 1/e al- 
titude is approximately 6000 feet, so tho density i3 approximately a factor 
of three smaller at 7000 feet than at 1000 feett 

Figures 3-9 and 3 — 1 0 show the time histories of contro's and tho ac- 
celerations at the pilot location. for 40 knots and 5 knots forward spaed at 
100 feet altitude. There is not much degradation in performance of the con- 
troller for 40 knot condition. There is, however, a slight degradation in 
performance of the controller at the 5 knot flight condition, but still the 
vibration is reduced to about 25? of the uncontrolled level. 

Variations in disc loading (helicopter weight) on the order of 20? did 
not significantly degrade the fixed-gain controller performance a 3 Indicated 
in Figures 3-11 and 3 - 1 2 . Though the steady- 3 tate controlled vibration level 
is roughly twice that of the controlled vibration at the design point, the 
performance i3 still acceptable. 
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Figure 3-K. Control Inputs and Accelerometer Outputs for Fixed-Gain 
Control at 120 kts/100 ft - Design Point 120 kts/100 ft 
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Figure 3-11. Control Inputs and Accelerometer Outputs for Fixed-gain 
Control at 120 kts/100 ft/AVI - +20$ (W-23520 lbs) 



Design Point 120 kts/100 ft 































3*3*3 Slnulatlon of Transient Flight Conditions 


The fixed-gain controller for the RSRA simulation manifested con- 
siderable robustness with respect to variation in flight condition and 
aircraft configuration (disc loading). In an attempt to demonstrate the 
robustness of the controller during non-steady-state flight conditions, the 
GENHEL simulation was set up to allow alteration of the trim control surface 
deflections. Since the SAS was not operative in the available GENHEL 
simulation program, the transient flight wa3 performed by interactive 
'manipulation of the stick', using ramp and step inputs for short periods of 
time while continually monitoring flight conditions. The command inputs 
used to accelerate from 80 knot3 to 125 knots forward speed and selected 
flight parameters are shown in Figure 3-13. Though the intent was to main- 
tain altitude at 1000 feet, this proved an exceedingly difficult task. The 
altitude time history indicates a loss of about 400 feet. The control feed- 
back commands and the vibrations for all three axes are presented in Figure 
3-1 4 . 

In these figures, no control is applied until 1 second. At 1 second 
the fixed-gain controller i3 turned on. From 3 seconds on, the input com- 
mand is applied. From 1.5 to 3 seconds the vibration level is reduced to 
less than 105 of its original level as is expected. During the transient 
flight, except during periods of radical command input, the vibration is 
reduced to about 15? of the uncontrolled level. Between 3 and 6.5 seconds, 
large control commands are being applied resulting in an increase in the RSS 
vibration level. Similar behavior is observed from 13 to 14 seconds during 
a period of Increased swa3h plate deflection During these periods however, 
the RSS vibration level does not exceed its unforced level. 
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Figure 3-13. Command Inputs and Selected Flight Parameters for Fixed-gain 
Control during Transient Flight from 80 kts/1000 ft to 
127 knot/595 ft 
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3. 3**1 Pl3CU33lon of Fixed-gain Controller Simulation Results 


Tatie 3-1 summarizes the results of the fixed-gain control applied to 
the GENHEL simulation of the RSRA. ha indicated earlier, their is a sig- 
nificant increase in the control authority required to achieve similar 
vibration reduction as the various flight conditions novo farther from the 
design point of 120 knots forward speed and 100 feet altitude. The fixed- 
gain controller performs quite well overall when applied to the simulation 
however. 


Table 3'1 . SUMMARY OF FIXED-GAIN CONTROLLER STEADY-STATE ACCELERATION AND 
CONTROL AMPLITUDES FOR STEADY-STATE FLIGHT CONDITIONS 


1 Flinilt Condition 

Steady-State 

Steady-State Amplitude Feedback 

Control 

VEQ 

Altitude 

RSS ACC 


(Peak-to-Peak) 


(knots) 

(ft) 

(ft/ 3CC 1 ) 

B l3 (deg) 

e Q (deg) 

A 1s (deg) 

O tr (deg) 






T" ] 

5 

100 

1.0 

0.3 

0.4 

0.4 

0.2 

40 

100 

0.9 

0.7 

0.4 

1 .0 

0.3 


100 

0.8 

1.0 

0.4 

0.3 

0.3 


1000 

0.9 

0.8 

0.4 

0.8 

0.3 

120 

2000 

1.0 

1 .0 

0.5 

0.5 

0.4 


5000 

1.5 

1.3 

0.9 

1.0 

0.6 


7000 

1.8 

2.0 

1.4 

1.5 

1.0 

80-127 

1000-600 

1 .5 

0.8 

0.4 

0.5 

0.3 




SECTION 4 


ADAPTIVE CONTROLLER FOR RSRA ROTORCRAFT VIBRATION SUPPRESSION 


This section discusses the issues involved in the design of an adaptive 
controller for the RSRA helicopter. The primary reason for considering an 
adaptive controller for vehicles of thi3 type 13 the wide range of input- 
output dynamic characteristics the helicopter is expected to assume as the 
flight conditions and/or structural configurations change. Rather than 
designing off-line, fixed-point control strategies, l.e. control strategies 
designed for a specific operating condition, and attempting to choose the 
appropriate strategy to meet the current conditions, the basic idea is to 
perform on-line system identification and control design simultaneously. As 
the operating conditions change, so does the system being identified and ap- 
propriate changes aro made in the control law. Thus, the control sy3tcm is 
adapting to the changing operating environment. 

Adaptive control sy3tem3 have two basic components: an 03timation and 
identification ’subsystem', and a feedback control ' subsystem' . The estima- 
tion and identification subsystem uses the inputs and outputs of the system 
in order to identify an equivalent input-output model of the system. In ad- 
dition, estimation of states in an internal description of the system to bo 
controlled may be performed in order to provide sufficient information to 
the controller subsystem. The distinction being made here is between full- 
state feedback versus output feedback control strategies. 

The controller subsystem uses the outputs from the estimation and iden- 
tification subsystem to compute feedback control commands. In a full-state 
feedback control design, an estimator i3 then used to obtain 'optimal' es- 
timates of the states in the dynamical model of the system based on past and 
present inputs and outputs of the system. Figure 4-1 gives a block diagram 
of the basic subsystems of an adaptive controller. The remainder of this 
section describes the design of the two basic subsystems and their 
Interconnections. Tho state estimation is performed utilizing a time-varying 
Kalman filter whoso system parameters are the outputs of a multi-input, 
multi-output parameter identification algorithm based on recursive predic- 
tion error methods. The controller is implemented via the time-varying LQR 
optimal control law which 13 dual to the Kalman filter estimator. 
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1M PARAMETER IDENTIFICATION AND STATE ESTIMATION 


State estimation and parameter identification is performed to provide 
to the control subsystem necessary information for the implementation of an 
effective control law. This is certainly not a simple task in the case of 
the RSRA helicopter duo to the extremely nonlinear nature of the vehicle. In 
this study, the state estimation and parameter identification were imple- 
mented sequentially rather than simultaneously. The option of augmenting the 
state estimation algorithm with states which are the parameters to be iden- 
tified was considered, but not implemented due primarily to computational 
considerations. This decoupling of state estimation and model parameter 
identification is 'suboptimal* however, and the possible performance 
degradation warrants further analysis. 

In thi3 subsection, the parameter identification algorithm which was 
employed i3 discussed in detail. The approximate model of the vibration 
dynamics is discussed, and the modeling and rejection of the periodic dis- 
turbances for identification purposes is detailed. Finally, discussions of 
practical issues such as implementation and tuning of the algorithm are 
given. 

4.1.1 Parameter Identification Algorithm 

The basic objective of parameter identification is to select fren a 
class of models, a model of the dynamics of the system under investigation 
which does the best Job, according to a given criterion, of predicting the 
future given the past. Implicit in the statement of the objective are three 
main tasks to be accomplished; selection of a class of models (model form), 
choosing appropriate criteria for discrimination between models in the 
class, and finally determination of the ’optimal' solution. The most dif- 
ficult or these tasks in general, i3 the selection of the model form. 
Models which are overly complex unnecessarily increase the computation load, 
and are more prone to suffering from input-output ldentif lability problems. 
On the other hand, models which are overly simplified may not retain suffi- 
cient information for adequate control. In designing the identification 
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subsystem, tho coal i3 to find the simplest model which results in an ade- 
quate oloaod-loop controller. Knowledge of the underlying physics of the 
problem, l.o. helicopter vibration, is used wherever possible. 

H. 1.1.1 Dynamic Model for Parameter Identification 

As discussed in the previous section, the majority of the power in the 
vibration is due to the excitation of the flexible modes of the RSRA 
helicopter by the periodic disturbances at n (N/rev) and its harmonics. A 
controller designed to suppress vibration at n and its harmonics will 
eliminate most of the vibration and could very well satisfy the objectives 
in terms of overall vibration reduction. 

The extremely narrowband nature of the vibration permits major reduc- 
tions in tho complexity of the model of the vibration dynamics. The 
majority of the vibration at the pilot location results from transmission of 
periodic disturbances at the rotor through the hub, amplified (or 
attenuated) and phase-shifted by the fuselage flexible modes. Therefore, 
models which can correctly account for the phase shift and attenuation frora 
each of the Inputs to each of the outputs in intervals of frequency around n 
and its harmonics should be sufficient. Linear second-order systeno satisfy 
these requirements for linear forced response. Thus, the model for ap- 
proximating the vibration dynamics amounts to a direct sum of underdamped 
second-order linear systems whose center frequencies arc the dominant vibra- 
tion frequencies, fl and its harmonics. A3 discussed in the previous 
sections, the inclusion of only the RSRA flexible fuselage modes near SI was 
intentionally designed to restrict the number of vibration frequencies to 
one (SI). 

Given a multi-input, multi-output system such a3 the RSRA helicopter, 
there are many model structures which incorporate the second-order single 
input to single output relationships desired. Noting that the flexible modes 
in the simulation were uncoupled by construction (cf. Section 2): and 
reasoning that for small vibrations the three orthogonal acceleration com- 
ponents would be nearly decoupled in an actual helicopter a3 well, three 
uncoupled multi-input single-output (MISO) models were used. With reference 
to equation 3«3, the equations for one of the axes can be written (in ob- 
server canonical 3tate-space form) as follows: 
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where x Is a two-vector [Xi , x s ]' of dynamic states, u is the vector of 
control inputs, r is the process noise/disturbance distribution matrix, and 
w is a vector of process nolse/disturb3nco terms: The subscript *p* is used 
to denote the 'plant. 

• It should bo noted that for performing vibration control, there are ac- 
tually four control degrees of freedom (inputs) assumed to be available, 
including the swashplato (two degrees of freedom), the main rotor collec- 
tive, and the tall collective. However, due primarily to it3 distance from 
the pilot location and its limited thrust capabilities in the lateral direc- 
tion, the tall rotor effectiveness in vibration reduction at the pilot 
location was found to be minimal. For the majority of the flight conditions 
investigated, the vertical, or z-comporiont of uncontrolled vibration ac- 
celeration wa3 the dominant component, and the lateral, or y-component was 
the smallest. Thus, in order to reduce the computation load and control ef- 
fort requirements, the tail rotor collective was subsequently ignored in the 
adaptive controller design. 



4. 1.1. 2 Disturbance Modeling and Rejection 


The vector w(k) in equation 4.1 is assumed to be composed of a broad- 
band, low power component modeling disturbances at frequencies other the £1 
and i ts harmonics (as well as accounting for dynamical model uncertainties 
to some extent), and a relatively high power narrowband disturbance at D. 
In order to obtain meaningful estimates of the parameters in the approximate 
second-order model of the Ir.put-output relationship, the effect of the nar- 
rowband disturbance terra which i3 present must be taken into account'. An 
important difference between the fixed-gain and adaptive control algorithms 
lies in the different approaches taken in addressing thi3 important issue. 
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In the off-line identification procedure U 3 ed in the fixed-gain con- 
troller do 3 ign, these disturbances were modeled as purely periodic with 
unknown amplitude and phase. The amplitude and phase were estimated by ex- 
plicitly including terms proportional to 0031^ and sirnj^ in the ARX models. 
Similar parameters could also be estimated in the on-line identification al- 
gorithm by employing model extensions similar to those used in the off-line 
identification procedure. A further extension would allow for estimation of 
narrowband disturbances by Including a damping coefficient in the model. 
This would seem appropriate in light of potential helicopter nonlinearities 
’broadening* the disturbance spectrum at a, The key point 13 that these ap- 
proaches require augmenting the helicopter vibrations dynamics model with a 
model of the narrowband disturbance process. 

The main objective in estimation of these periodic disturbance terms, 
however, is basically to eliminate the disturbance (’colored nol 3 e* ) com- 
ponent from the measurements. Thl 3 i 3 done prior to updating the parameters 
in the identification algorithm, otherwise biased parameter estimates result 
(leading to potential problems in the controller). Realizing the objective 
is to eliminate a narrowband disturbance of known center frequency, a com- 
putationally more attractive, but 3 uboptimal approach was taken. The 
elimination of the frequency components at N/rev was accomplished with a 
digital notch filter. The continuous domain transfer function equivalent of 
the digital filter employed ha 3 two zeroes on the imaginary axis at N/rev, 
and two poles Just inside the left half-plane: 


F BP (3) 


2 2 

s + n 

2 2 
a + 2efls ♦ n 


As a practical issue, in order to improve the transient performance of 
the notch filter for parameter identification and retain the narrow 
bandwidth required during the identification process, a time-varying damping 
coefficient was implemented. During the transient region of the simulation 
start-up, a relatively wide bandwidth notch (c~2l)’ is used. The damping 
coefficient is decreased exponentially when parameter identification is in- 
itiated, narrowing the notch filter bandwidth significantly (c- 0 . 2 ?). Note 
both the inputs and outputs are notch filtered before they are passed to the 
parameter identification algorithm so as not to introduce an distortion of 
one relative to the other. 
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To conclude the discussion of the model of the vibration dynamics and 
the disturbance modeling (rejection), several important points concerning 
the number of identifiable parameters (to be estimated) in the model in 
equation 4.1 are noted. Without loss of generality, the 3tate component of 
the output was chosen to be the first state x, (i.e. observer canonical 
form). Since a second-order system is completely specified by its natural 
frequency and damping coefficient, only two parameters are required in the 
system dynamics matrix. The remaining elements in the control distribution 
matrix are unconstrained. By notch filtering the measurements and inputs, 
identification of parameters in the disturbance distribution matrix as well 
as possible parameters in a dynamical model of the narrowband disturbances 
is not required. 

4. 1.1. 3 MISO RPEM Parameter Identification Algorithm 

Decoupling of the model into three separate axes led to considerable 
simplifications in the design of the parameter identification scheme. The 
MISO state-space model in equation 4.1 i3 readily transformed into an equi- 
valent MISO ARX (autoregressive model with exogenous Inputs) model, which is 
a standard form used in several recursive prediction error parameter iden- 
tification (RPEM) algorithms. Transforming equation 4.1 into an ARX form 
gives: 

y(k+1) - A 1 y(k)+A 2 y(k-1)+B Q u(k+1 )+B 1 u(k)+B 2 u(k-1 ) + v(k) (4.2) 


where in an obvious notation, the arguments are time indices and y is a 
scalar output with u the vector of control inputs. The coefficients are 
matrices of the appropriate dimensions, and bear the following simple 
relationships to the parameters in the state space model: 


F 

F 

G 


11 

21 
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B ij * A i B 0J 

• B or 


(4.3) 
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By stacking the outputs and inputs on the right-hand side of equation 4.2, 
and collecting all the parameters in a single vector, the equation can be. 
writton in regression form as indicated in equation 4.4: 

y(k) - 0 T <»(k) + v(k) . (4.4) 

A derivation of the classical recursive least-squares Identification 
algorithm is outlined below, leaving detailed derivations to the references. 
The least-squares cost criterion on the parameter vector estimate is: 

N T ? 

J “ Jo Cy(k) ' 0 ♦<.*>] ■ ■ 

Minimization of this cost function is accomplished by differentiating with 
respect to the parameter vector (0) and setting the gradient to zero. 
Taking care to preserve the recursive nature of the equations, the recursive 
least-squares identification algorithm is obtained. 

R(k) - Fl(k-I) ♦ -i-CMk),> T (k) - R (k— 1 ) ] 

0(k) - 0(k-1) ♦ ~ R(k)"V(k)c(k), (4.5) 

e(k) - y(k) - 0 T (k-1)$(k). 


Though these equations are recursive, they are not well suited for on-line 
implementation 3ince at each 3tep a matrix inverse is required. However, 
taking advantage of the fact that the update to the information matrix (R) 
is of rank one, the matrix Inversion lemma can be used to exchange the 
matrix Inverse for a scalar division in the propagation of the parameter 
covariance matrix as shown in equation 4.6. 

P(k) - ?(k~1 ) - (k ~ ' (4.6) 

1 + <f> (k)P(k-1)<>(k) 
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The parameter update equation becomes: 

K(k) - P(k-1)<J>(k)[1 + <> T (k)P(k“1 )$(k)] 
9(k) - 0(k-1) ♦ K(k)e(k) . 


(4.7) 


The above equations assume that the parameters are constant over time, 
though the estimate e(k) is certainly a time-varying quantity. The inter- 
pretation to be given is that 0 ( k > is the best estimate of the conr'.unt 
parameter vector given all the past information. In order to allow for 
time-variation of the parameters, past information must somehow be 
deweighted and greater emphasis placed on more recent information. 

There are several methods for modifying the above equations to accom- 
plish thi3 goal. The first, and computationally less intensive method is to 
insert a ' forgetting, factor' into the cost function such that more weight is 
placed on recent errors than past errors. If the weighting is exponential, 
then above equations are simply modified: 


POO - JL PC*-!) - 

K } X(k) * <>(k)P(k“1 )<? (k) 

(4.8) 

K(k) - P(k-1)<>(k)/(X(k) + <j>(k)P(k-1)4> T (k)), 

where X(k) < 1 is the 'forgetting factor', a scalar variable (usually set to 
some fixed value). These equations effectively Include a time update along 
with the measurement update by scaling up the parameter covariance matrix P 
by the factor 1/X(k). 

There is a certain amount of redundancy in equations 4.8, and this fact 
can be taken advantage of in designing a fast algorithm for computing the 
necessary quantities. The fast implementations rely on the fact that a 
positive definite matrix (all covariance matrices are positive definite) can 
be factored into a product of three matrices; an upper (lower) triangular 
matrix with unit diagonal elements, a diagonal matrix w ith positive ele- 
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dents, and the transpose of the upper (lower) triangular matrix, the so- 
called' U-D algorithm. These equations propagate an upper triangular factor 
(U) and a diagonal matrix (D) such that 

P - UDU T . 

The U-D update equations are basically a modified version of the Gram- 
Schraidt orthogonalization procedure, the details of which are left to the 
references (cf. Bierman, [9]). 

Returning to the problem of the time update of the parameter co- 
variance, the fact that only a scalar forgetting factor is allowed i3 quite 
restrictive in the sense that all parameters are assumed in some sense to 
have the same time-varying time constant. In most practical situations, ex- 
cellent arguments based on physical laws can be made to suggest that some 
parameters have larger expected time-variation than others. An identifica- 
tion algorithm which does not take this information into account must cer- 
tainly suffer performance degradation in comparison to an algorithm which 
does. 

The discussion in the preceding paragraph suggests that an implementa- 
tion of the covariance time update which allows for the inclusion of more 
(stochastic) information is more appropriate. This second method of 
deweighting past information is more appropriately thought of as propagation 
of the filtered parameter covariance in the presence of additive disturb- 
ances (process noise or model uncertainty). The time update equation is now 
separate from the measurement update and is given by: 

P(k+l|k) - P (k | k) + Q(k) (4.9) 

where the notation has changed to reflect the fact that at each time k there 
are predicted and filtered covariances. Though the computation load i3 in- 
creased, there is a modified weighted Gram-Schnidt procedure (cf. Bierman, 
[9]) for performing this time update in the U-D factorization formalism. 
For the reasons already cited, this method of covariance time update was 
chosen over the 'forgetting factor' approach. Tests of both approaches in- 
dicated that the choice was well-founded on the basis of robustness and 
overall vibration suppression performance. 
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ll. 1.1.1* Practical Considerations for ID Algorithm Implementation 1 


An important consideration in parameter identification i3 that of Input I 

design. Identification of parameters in a closed-loop system without per- ! 

sistently exciting external inputs is a singular problem in the sense that 
parameters in the open loop transfer function are not identifiable. In or- 
der to insure a well-posed identification problem, a sufficiently wide 
bandwidth external input is required to excite all the modes of the system. 

i 

For identification of the parameters in the narrowband model of the RSRA j 

helicopter dynamics near N/rev, such an input was constructed by passing j 

random noise through a cascade of two second-order filters whose damping ! 

coefficients were 30? and whose center frequencies were fl. Figure i)-2 gives 
a two-second sample of this excitation as well as a plot of it3 power 
spectrum. 
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Figure l»-2. Sample Time History and Power Spectrum of Wideband 
Excitation for Identification. 


The issue of how much noise to add for identification basically in- 
volves a trade-off between the requirement for persistent excitation at a 
signal-to-nolse ratio (SNR) sufficient to provide reliable parameter esti- 
mates and the desire to keep the control effort and. induced vibration out- 
side the controller bandwidth below acceptable levels. In the case of the 
RSRA helicopter, acceptable control surface deflections for vibration con- 
trol (at Il/rev) were assumed to be on the order of two degrees for all the 
controls. 

In the implementation, the noise amplitude was made a function of the 
performance of the control) er. The controller performance measure used was a 
sum of low-pas3 filter outputs whose inputs were the rectified accelerometer 
outputs from each axis. Between preset upper and lower limits, the input 
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noise amplitude .was a linear function of the difference between the perfor- 
mance measure and a preset ID-OFF threshold (discussed in the next 
paragraph). This had tho effect of increasing the input amplitude when the 
performance measure indicated a need for improved parameter identification, 
and reducing the control deflection required to suppress vibration resulting 
from this noise when the identified parameters were sufficient for control 
purposes. 

The performance measure was used not only to regulate the wideband ex- 
ternal input amplitude, but also to turn on and off the identification 
algorithm. Using hysteresis logic, ID was turned off as the performance 
measure fell below a preset threshold (3 ft/sec 1 is used for most of the 
results presented in this section), and was turned on if the performance 
measure exceeded a second preset threshold which was larger than the first 
(3.6 ft/sec 1 was used for most of the results presented in this section). 

This switching logic is primarily used to prevent the identification 
algorithm from diverging during periods when identifiability of the’ 
parameters is low. Parameter identifiability decreases as tho amplitude of 
the external wide bandwidth noise input decreases ,. and the parameter iden- 
tification algorithm becomes singular when- the exogenous wide bandwidth 
input is removed. Since the controller bandwidth is narrow, it is desirable 
to eliminate the wide bandwidth noise when a steady-state is reached. This 
enabling/ disabling logic prevents the burst-type of instabilities reported 
in tho literature for adaptive controllers. 


4.1.2 State Estimation Algorithm 

The state estimation algorithm uses the model of the system to be con- 
trolled which, as discussed In the previous section, is comprised of three 
decoupled second-order systems. The parameters in these second-order models 
are functions of the outputs of the parameter identification subsystem, and 
are given by the relationships in equation 4.3. The states of tnese second- 
order models are estimated using the measurements, and these estimates are 
used in the feedback control law. ' In the limit of zero measurement noise on 
the accelerometer outputs and control surface deflections, state estimation 


- 60 - 



reduces to bnndp333 filtering of the outputs aa was employed In the fixed- 
gain algorithm. The inclusion of a state estimator in the adaptive 
algorithm makes the feedback controller less susceptible to input process 
noise, input measurement noise and output measurement noise at the C03t of 
an increase in the computational load. 

4. 1.2.1 Dynamic Model/Frequency Shaping Filters for State Estimation 

The second-order model which approximates the dynamics of the RSRA 
hollcoptor around N/rev has been discussed in detail in the previous sec- 
tions. The state estimation algorithm augments thi3 model with the 
frequency shaping filter states as discussed in detail in Section 3.2.1. 
The parameters associated with these artificial states are calculated based 
on the frequency shaped filter requirements and are not parameters to be 
identified. These frequency shaping filters effectively penalize the 
bandwidth of interest in the cost functional in the LQR problem formulation, 
and can bo thought of as simply providing dynamically compensated outputs 
for feedback control. Only that frequency v/hich is dominant in tho vibra- 
tion spectrum is fed back through the controller so a3 not to excite any 
other vibration modes via helicopter nonlinear! ties. 

Tho dominant frequency in the RSRA helicopter (simulation) vibration 
spectrum is ll/rev and is nearly constant. Variations in angular velocity 
are expected only under tavere, or rapid, control demands, possible 
regulator failure (a rotor speed regulator was included in the simulation), 
and random disturbances. The random disturbances are expected to bo small, 
but important in the sense that the phase of the disturbance terms is 
altered. These observations formed the guidelines for the design of the 
frequency shaping filters. 

The frequency shaping filters used in the control design have conflict- 
ing design requirements. Lesa control effort will be used in suppressing the 
H/rov vibration a3 the filter bandwidth is decreased. However, Infinite Q 
(or zero bandwidth) filters are not desirable from the standpoint of either 
robustness or speed of response of the adaptive system. On the other hand', 
relatively large bandwidths (e.g. 5? of N/rev) excite modes which are out- 
side the region of validity of the second-order model used, and consequently 
may not be controllable. Values of the damping coefficient from 0.1? to 0.5? 



VST" • 


gave satisfactory response times while sufficiently suppressing frequencies 
other than N/rov. 

In an attempt to employ a strategy similar to that used in the notch 
filter design (cf. Section 4.1.1. 2), the capability to vary the damping 
coefficients as a function of time and controller performance was also 
implemented. A wider bandwidth i3 used during the transient region to im- 
prove the response time, followed by a period of decreasing bandwidth to 
improve the steady-state vibration suppression performance. Damping coeffi- 
cients on the order of 0.5J were used during the transient region, 
exponentially decreasing to 0.052 after sufficient vibration suppression was 
achieved. This resulted in improvement in the 3teady-state performance, 
however it was not significant enough to warrant the increased computational 
load (improvement was on the order of 3dB over fixed damping coefficient 
filters). 

Several other filter design issues were also addressed. Zero phase 
shift at W/rev effectively eliminated any estimator lags that the controller 
would have to compensate for with leads (differentiation). Care was taken 
in tho digital mechanization to accurately determine the center frequency of 
the filters. Finally, normalization of the filter gain at tho center fre- 
quency reduced round-off error propagation in tho estimator and controller 
equations. 

Based on the frequency shaping filter design requirements, a filter 
with a zero at the origin, unity gain at Q , and complex pair of poles deter- 
mined by the desired damping coefficient and a natural frequency of n was 
constructed. Assuming a measurement of the rotor angular velocity was 
available, the filters were continuously tuned to a center frequency of 
N/rev, tracking any variations in the rotor frequency. The bilinear trans- 
formation from continuous to discrete domain with prewarping was used, 
giving tho desired digital filter. The critical frequency for prewarping 
was n. 

In the continuous-time domain, the frequency shaping filters have the 
following Laplaca transfer function; 

2chs 


H(S) 


s 2 + 2c0s + n 2 
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The transformation to dl 3 crctc-time domain first involves substitution of n 
with £1^, the prewarped critical frequency: 

. 2 -AT. 

tip - - tan(-) 

Then the Laplace variable 3 is replaced with; 

. 2(z - 1) 

m * i) 

where T 13 the sampling interval. Expanding the resulting numerator and 
denominator expressions results in the following discrete equivalent trans- 
fer function: 


H(z) 


V * V * b 2 
V 2 + V + a 2 
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Note that this transformation results in a filter with direct feedthrough 
since the numerator and denominator polynomials are of the same order. This 
results in a requirement to modify the standard control gain calculations as 
was discussed in Section 3.2. Realizing this transfer function in observer 
canonical state-space form gives the following 3econd-order system for the 
frequency shaped filters: 


x(k«-1) - 
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fs21 
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x(k ) + 


fall 


f 321 


y(k) , 


z(k) - [1 0] x(k) + D fs y(k), 
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The output of the frequency shaped filters (z) is assumed with loss of 
generality to be the first state plus a direct feedthrough term. However, 
since the controller feeds back the states of the model only, the filtor 
output need never bo constructed. The effect of the direct feedthrough term 
is to modify the costs in the performance function for the controller design 
(cf. Section 3.2 and Appendix C). 

4. 1.2. 2 Measurement Models for State Estimation 

The accelerometer outputs can be modeled without I033 of generality as 
the first of the two states in the second-order models for the input-output 
relationships in each axi3. A direct feedthrough term is also estimated in 
the ARX form of the transfer function, but since significant cost is placed 
only on the frequency shaped filter 3tates, the sole effect of the direct 
feedthrough torn is to modify the elements- in the augmented control distri- 
bution matrix G. 



Figure *4-3. Frequency Response of Bias Rejection Filter 


The issue of biases in the accelerometer outputs must be addressed. 
Biases can either be estimated (identified along with the other parameters 
in the model) 30 there effect can be removed, or more directly, removed from 
the measurements before entering the estimation and identification sub- 
system. Since the frequency range of -interest for vibration control is well 
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above DC, a first-order lead-lag filter with a zero at DC and a pole at 1 Hz 
is sufficient. Transforming this filter into the discrete domain U3ing the 
zero-order hold equivalence transformation gives the following discrete 
transfer functions 

F(z) “ z “ 0.9752 * 

assuming a sample rate of 250 Hz. A frequency response plot of the bias 
rejection filter is shown in Figure *4-3. 

4. 1.2. 3 Kalman Filtering for State Estimation 

Having a set of three decoupled fourth-order systems (a second-order 
approximate dynamics model with a second-order frequency shaping filter) 
with parameters predetermined as well as identified in the parameter identi- 
fication subsystem, state estimation is performed using a Kalman filter. _ 
The dynamical model of one of the fourth-order subsystems is given by: 



r* ** 

F 0 
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x(k+1) - 

> H p F fs. 

x(k) + 

G. D 
fs p_ 


u(k) ♦ w(k) 


y(k) - [1 0 0 0]x(k) + DpU(k) + v(k) , 


where the elements of the system matrices have been previously defined. 
Denoting by F the entire 12 by 12 system matrix, and by G the entire 12 by 3 
control distribution matrix, the Kalman filter equations can be written: 


x(k+l|k) - Fx(k[k) + Gu(k) , 
v(k+t) - y(k+1) - Hx(k+l|k) , 
x(k+l|k+1) - x(k+l|k) + K(k+1)v(k+1) , 
P(k + 1 | k) - FP(k |k)r T + Q(k) , 
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R v (k>1) - HP(k+1 |k)H T * R , 

K(k+1 ) - P(k+l|k)H T R v (k+1)" 1 , 

P(k+1 |k+1 ) -[I - K(k+1 )H]P(k+1 |k) , 

where H - [1 0 0 0] is the measurement distribution matrix, R is the 

covariance of the measurement noise (v), and Q is the covariance matrix of 
the process noise (w). The process noise is included in the estimator to 
model the source of the periodic disturbances (the process noise is filtered 
through the bandpass filters) as well as to account for the approximate na- 
ture of the dynamical model. P is the covariance matrix of the state 
estimation error, and K is the Kalman gain matrix which corrects the state 
estimates (x) based on the prediction error, or innovations (v). 

To conclude this subsection on state estimation, a few algorithm im- 
plementation issues are addressed. The Kalman filter equations ' were 
implemented in square-root form to decrease storage requirements and com- 
putational load, as well as to ensure the positive definiteness of the 
covariance matrix. The details of the 3quare-root formulation can be found 
in the references (cf. Bierman, [ 93 ). 

The filtered state estimates are the estimates actually passed to the 
controller for feedback control computation. If this estimator/controller 
combination is mechanized in a real-time environment, the control law compu- 
tation by a flight computer will result in a delay in the control actuation. 
Further study to determine the effects of the delay on the stability of the 
controller will be required. 

4. 2 CONTROLLER SUBSYSTEM 

U .2.'. Adaptive Controller Gain Calculation 

The design of the controller subsystem utilizes basically the same 
linear quadratic regulator (LQH) design procedure discussed in detail in the 
previous section. The main difference is that the dynamical model of the 
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system i3 tino-varying, thus requiring on-line gain computation. Also as 
discussed in the previous section, the standard LQR problem formulation must 
be modified to include the effects of a component of the Inputs in the 
outputs. This modification i3 described in detail in Appendix C. 

Using the dual formulation of the estimation problem discussed in the 
previous subsection, the same equations in square-root form can be used 
(with appropriate variable substitutions) to obtain the optimal time-varying 
controller gains. Basically, the measurement distribution matrix is re- 
placed with the control distribution matrix, and the process and measurement 
noise covariance matrices are replaced with state and control cost matrices 
respectively. The appropriate substitutions can bo found in the references 
(cf. Bryson and Ho, [8], Franklin and Powell, [9]). Using the time-varying 
'optimal' gains gives the following control feedback law: 

u(k) - C(k)x(k|k) , 

where C(k) i3 the optimal Ricatti controller gain, the dual of the Kalman 
filter estimator gain K(k), suitably modified to account for the direct 
feedthrough term a3 discussed in Appendix C. 


k.2.2 Determination of State and Control Cost Matrices 

The parameters which must be predetermined for control gain calculation 
other than the frequoricy shaping filter parameters are the state and control 
cost matrices. Following- the heuristic rules of thumb proposed by Bryson 
and Ho, the control costs are set to the inverse of the 3quare of the maxi- 
mam RMS control authority desired. 

I* 5 peak-to-peak deflection »> 0.707° RMS -> cost of 2 (deg 2 ) 

The state costs are similarly set to the inverse of the square or the 
maximum steady-state RMS vibration level acceptable. Cost is placed only on 
the frequency shaping filter output states to insure that the majority of 
the feedback control authority is at fl. 
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In order to provent transients in the controller, the state costs are 
time-varying as well. When the controller is first turned on, the initial 
Hicatti matrix (tho dual to the state estimate error covariance matrix) is 
set nearly to zero and the state cost3 are set to zero. The costs are then 
exponentially increased to their final predetermined level. The exponential 
time constant was iteratively determined to give satisfactory transient con- 
troller performance, the primary objective being prevention of large control 
deflections from exciting tho nonlinearities of the RSRA helicopter 
simulation. Since the Ricatti equation 13 nonlinear, the transient perfor- 

mance is a function of the magnitude of the state costs as well as their 
time- dependence. 

To summarize this discussion of the adaptive control design, a more 
detailed block diagram of the controller is given in Figure 4-4. 

4.3 ADAPTIVE CONTROLLER TUNING AND SIMULATION RESULTS 


The adaptive controller discussed in detail in this section was pro- 
grammed and added to the GENHEL simulation of the RSRA helicopter. In thi3 
subsection, he results of simulations run at various operating conditions 
are presented. The initial conditions for the runs at each operating condi- 
tion are the same as 13 the tuning of the .'ariou3 process noise covariance 
matrices and regulator cost matrices. Since these matrices represent in a 
stochastic sense the uncertainty in the dynamical model and relative costs 
of control and acceleration at N/rev respectively, they could easily be 
scheduled on dynamic pressure, rotor angular velocity, altitude, or 3ome 
other appropriate parameter if further study showed that there is a corres- 
ponding change In the model uncertainty or control and state cost3. 
However, no attempt to optimize the performance at each operating condition 
over tho various weighting matrices (process noise and cost matrices) was 
made in the analysis presented here. 



















>J . 3 . 1 Initial Conditions and Algorithm Tuning 


The Initial conditions were the 3ame for each of the run3 at the 
various operating conditions. The Initial parameters for the three 
decoupled oscillators which modeled the dynamics of the helicopter around 
N/rev were the same: 

A, ■ 1.7 , 

A 2 - -0.9 . 

1 6 — -j 0 

Recall for a pair of discrete poles at pe J and pe J , 

A. - 2pcose , 

a 2 - -p 2 . 

Thus, these initial conditions represent stable poles (inside the unit 
circle) noar n (at a sampling rate of 250 Hz) and lightly damped. The ini- 
tial values of the entries in the control distribution matrices (B) are all 
zero. The initial covariance of these parameters Is diagonal with a value 
of 10 for each parameter variance In the parameter identification algorithm. 
Moderate changes in the initial covariance do not affect the final results 
significantly for two reasons. First, after allowing the GENHEL simulation 
internal states one second to reach steady state, the identification algo- 
rithm i3 given another second to accumulate Information before the 
controller is turned on. Secondly, process noise is included In the model 
which deweights the contributions of initial conditions as data are 
processed. 

The initial values of all the states, in the estimator are zero, and 
their corresponding covariance matrix is the identity matrix. Since the es- 
timator is turned on even before the identification Is enabled, and since 
process noise is included in the model for these states, large changes in 
the initial state covariance will net affect the final results signifi- 
cantly. 

The process noise variances on the parameters ana the states in the es- 
timator were determined iteratively. At >40 knots and 1000 feet altitude, 
these variances were iteratively changed until satisfactory performance was 
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achieved. (As an aside, only two iterations were required, since excellent 
initial values could be calculated based on expected parameter variations 
over tine. A priori information such as this can be used to schedule these 
variances in the future if desired.) The process noise variance associated 
with the parameters was required to allow, the parameters to vary with time. 

The equivalent dynamical model for the parameters is: 

• 

e - w 

where w is white Gaussian noise of variance density corresponding to the 
discrete excitation noise variance (a factor of the sampling interval is 
Involved). For the oscillator parameters (Aj.Ag), the variance used wa3 

0.001, and for the remaining parameters (B), the variance was 0.01. The 
variance of the oscillator parameters is unitless, while that of the para- 
meters in the control distribution matrices is in units of degrees squared. 

(The continuous variance density would have units degrees squared per 
second.) 

The process noise variance used for the estimator states was the iden- 
tity matrix. This value, though somewhat arbitrary, is reasonable based on 
the fact that the states in the estimator are in effect normalized 30 that 
their magnitudes remain near unity. Though the magnitude of the process 
noise variance is important in determining the covariance of the state es- 
timation error, its value relative to the measurement noi3e variance is the 
key factor in determining the optimal Kalman estimator rains. A value of j 

0.01 ft/sec 2 for the standard deviation of the accelerometer noise wa3 used 
for the results presented in this report in which measurement and process 
noise were absent. In cases where measurement noise was added, an ap- 
propriate measurement noi3e 3igma wa3 used in the estimation algorithm. j 

The cost matrices for the regulator gain calculations were also deter- j 

mined iteratively. Initial values were calculated by the Bryson method I 

mentioned earlier, and fine tuning to obtain acceptable performance was per- j 

formed at MO knots and 1000 feet altitude. The final values for the 3tate t 

and control cost3 used for the results presented herein are 5 for the states I 

representing the frequency shaping filter outputs, zero for the remaining j 

states, and a cost of 2 on all three controls; A, , B, , and 8 . I 

ts Is c | 


« 




As mentioned earlier, the controller i3 turned on after two seconds to 
allow transients in the GEfJKEL simulation and the estimation and identifica- 
tion algorithms to settle. The rate at which the control gains approach 
their steady-state values (assuming for the moment that the parameter varia- 
tions are insignificant) is determined by two factors; the initial value of 
the Rlcatti matrix (the control dual to the state estimate error covariance 
matrix), and the time history of the state and control cost matrices. The 
Rlcatti matrix is initialized to an extremely small factor times the iden- 
tity matrix. This results in essentially zero initial control gains (yet 
avoids numerical difficulties which would arise if zero were input 
instead!), a desirable feature in preventing po33ible transients. The expo- 
nential growth of the state costs from zero to their final fixed value also 
helps prevent undesirable transients in the control demand. The time con- 
stant of this exponential growth is approximately two seconds for the re- 
sults presented in this report. 

It should be emphasized that these intial parameter values are required 
for the most part to overcome the phenomenon of algorithm start-up , not 
helicopter start-up . Once the RSRA helicopter and the adaptive controller 
are past. the initial transients, if the algorithm is properly tuned, it will 
track the system variations and should not require reinitialization. How- 
ever, should it become necessary to reinitialize the algorithm as the result 
of a detected Instability, the same initialization procedure is engaged. 

As discussed in the previous subsection, once the algorithm has 
achieved an acceptable 3teady-state performance as determined by a cumula- 
tive acceleration magnitude threshold, the parameter identification 
algorithm is disabled in order to alleviate the need to add noise to the 
control Inputs (for identifiabllity) and to reduce the computation load. 
The performance Is continually monitored, however, and should it degrade 
past a second cumulative acceleration threshold, the identification algo- 
rithm is re-enabled. When the identification algorithm is re-engaged, the 
parameter covariance is augmented by adding 0.1 to all the diagonal 
elements. Since the models of parameter dynamics include process noise, 
this increase in the covariance serves merely to speed up the adaptation 
time constant. The overall tracking performance is still determined by the 
process noise variance, not by the initial parameter covariance. 



U . 3 . 2 Adaptive Controller Simulation Results 

This subsection presents graphically the results of the adaptively con- 
trolled GSNHEL simulation -uns of the RSRA helicopter at various flight 
conditions. The plots of accelerometer outputs are the outputs of the 
GENHEL simulation program, filtered through a high-pass filter as discussed 
In section *4.1 ,2.2 to remove any low frequency components. Removal of most 
of the rigid body accelerations yields plot3 which are conveniently nearly 
zero-mean. Also, unless otherwise indicated, the plots of control surface 
deflections are plots of the feedback correction components only, and do not 
include any non-zero set point commands from a stability augmentation sys . 
(SAS) or the trim condlt-ons. Thus, the plots of control deflections are 
also nearly zero-mean. 

The time histories plotted clearly manifest the controller time con- 
stants as envelopes of the H/rev modulation and its harmonics. Also, aa 
discussed previously, the first one second of output is uusuppre3sed vibra- 
tion of the RSRA helicopter. Thus, the first second of the accelerometer 
outputs gives an indication of the unforced vibration levels (at the pilot 
location) for each flight condition. Furthermore, during the second one- 
second Interval of the output time histories, the effect of the additive 
noise can be seen by comparing the vib-at.ion levels during the first and 
second one-second intervals. The control del lections during the second 
one-second Interval. Indicate the amplitude of the additive excitation 
clearly. 

The RSS acceleration output presented in each f igure is the root-sum- 
square of the three orthogonal accelerometer outputs shown. Low frequency 
rigid body acceleration ha3 been filtered out, so the RSS acceleration is a 
direct measure of the total vibra' ion power at the pilot location. 

Figures 14-5 through *4-9 give selected results of t ,- e adaptive con- 
troller for velocities of 5 knots, HO knots, and 1 ?n knots at an altitude of 
100 feet. Plots of the time histories of selected parameters and their as- 
sociated estimate error sigmas for the *40 knot/100 feet altitude flight 
condition are included as well. The notation used for the labeling of the 
ordinates in the parameter time histories indicates the input-output pair 
(eg. AZ/B1S) in the case of control distribution parameters (B0.B1), and the 
output only (eg. A/) in the ca3e of the dynamic model parameters (.k' ,k2) . 
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Table 4-1 gives the correlation matrices for the three decoupled MISO 
parameter covariance matrices for’ the 40 knot3/100 feet flight condition. 
The correlation between the f^and A,, terns in each axis is approximately - 
0.9 as is the correlation between the B Q and B 1 terms in each of the input- 
output pairs (cf. [3,6], [4,7], and [5,8] elements of the correlation 
matrices). This is a consequence of the strong identif lability of the 
input-output amplitude ratio, and the comparatively weak identif lability of 
the relative phase at Cl. 

Figures 4-10 through 4-12 present selected results of the adaptive con- 
troller at 5 knots, 40 knots, and 120 knots at 1000 feet. Power spectra of 
the controlled and uncontrolled steady-state vibrations in each axis at 120 
knots and 1000 feet altitude are shown in Figure 4-13. These spectra indi- 
cate the the vibration at Cl has been suppressed to levels below the 
contribution of the higher harmonics and indicatos the adaptlvo controller 
is performing as expected. Figures 4-14 through 4-17 present the results of 
tho adaptive controller for a speed of 120 knoto and altitudos of 2000, 5000 
and 7000 feet. 

The results in Tlgures 4-5 through 4-17 indicate that from the same al- 
gorithm initial conditions, the adaptive controller converges at each flight 
condition, yielding excellent performance lr» terms of RSS vibration suppres- 
sion and control effort expended. To investigate the performance of the 
adaptive controller under changing flight conditions, two cases were run. 
Figures 4— 1 8 through 4-20 show the results for a case in which the Cl (N/rov) 
was changed abruptly by decreasing it 10 rad/sec in 0.2 seconds, held con- 
stant for the next 1.8 seconds, and then returned to its original value 
linearly over the next 4 seconds. The objective of this test was to inves- 
tigate the ability of the adaptive algorithm to track variations in rotor 
angular velocity and resulting changes in the helicopter dynamics. 

The objective of the second case was to investigate the performance of 
the adaptive algorithm during a transition from one flight condition to 
another. The initial flight condition U3ed for this test was 80 knots at 
1000 feet. Control commands were given to Increase the velocity to 125 
knots. Figures 4-21 through 4-23 present selected results for this case, 
including plots of trajectory parameters. 
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Table it-1 . IDENTIFIED PARAMETER CORRELATION MATRICES FOR 
ADAPTIVE CONTROL AT AO KHOTS/IOO FT 


X-axlo Parameter Correlation Matrix 









A1 

A2 

BIS/BO 

A1S/BO 

COL/BO 

B1S/B1 

A1S/B1 

COL/BI 



i • GC5D 

-o.tm 

-0.2172 

-0.4523 

0.1034 

0.6352 

0.S123 

-0.2SE0 



-0.0213 

urn 

0.2731 

0.3107 

0.0232 

-0.2312 

-A 

0.1034 



-0.2122 

0.2731 

1.C3G0 

0.2044 

-0.2574 

-0.0353 

-0.K33 

0.0127 



-0.0323 

0.2107 

0.2044 

1.0389 

-o.iesi 

-0.O423 

-0.0344 

0.1427 



0.20 

0.C2S2 

-4.2S74 

-0.1631 

1.0350 

0.45S) 

0.2521 

-4.033? 



0.CSJ2 

-0.£-3i2 

-4.K31 

-0.0423 

0.4553 

1.6553 

0.1089 

-0.2377 



0.GI23 

-0.S233 

-0.3503 

-0.C3J4 

0.2321 

0.1043 

i.C3» 

-0.2233 



-0.2053 

0.1034 

0.01S7 

0.K57 

--O.C337 

-0.2377 

-0.2253 

1.0223 





Y-axla Parameter 

Correlation Matrl 

X 



I 

A1 

A2 

BIS/BO 

A1S/B0 

COL/BO 

B1S/B1 

A1S/B1 

COL/BI 



1.03 

-0.C373 

-0.G1C5 

0.1727 

-S.7D23 

-4.0281 

0.0723 

0.1610 



-0.C370 

i.c-::3 

o.cr.o 

-0.2553 

6.1353 

-Mm 

0.1037 

-4.002-3 



-O.CICj 

0.0533 

l.CSS) 

0.1707 

-0.2937 

-0.&S9 

-0.MJ 

-0.0D4 



0.1227 

-S.2365 

0.1737 

1.CSS9 

-0.1153 

-0.6352 

-0.0731 

0.0223 



-0.2323 

0.2323 

-0.2237 

-0.1153 

1.0GC9 

0.KS3 

o.iss 

-0.0327 



-4.CS31 

-O.C4G4 

-o.Gsr 

-0.0252 

o.5f:j 

1.6059 

0.1702 

-0.2723 



0.0253 

0.1437 

-0.2533 

-0.0731 

0.1825 

0.1702 

1.6223 

-0.1052 


t 

l 

0.1414 

-3.2043 

-0.0254 

0.0313 

-0.6327 

-0.2723 

-Q.m> 

1.C899 


J* 

i 

{ 



Z-axla Parameter 

Correlation Matrix 



» 

1 

i 

A1 

A2 

BIS/30 

A1S/B0 

C0L/30 

B1S/B1 

A1S/B1 

COL/BI 


| 

i.mo 

-0.0722 

-0.2773 

0. CS23 

-0.1303 

0.2089 

-0.2072 

-4.CE27 


[ 

1 

-0.0722 

1.C003 

0.222:,' 

0.0460 

0.2525 

-0.2332 

0.0223 

-0.1102 


i 

i 

-0.277S 

0.2223 

l.CMO 

0.1710 

-0.2052 

-0.0253 

-9.1033 

0.0334 


} 

0.0323 

0.0403 

0.17IC 

l.CWO 

0.0553 

-0.1457 

-0.0753 

-0.1G10 


1 

| 

-0.1303 

0.2K5 

-0.2532 

0.0SS 

1.6000 

0.4811 

0.0243 

-0.CM2 


1 

> 

0.2430 

-0.2522 

-0.K33 

-0.1407 

0.4011 

1.6550 

C.18S7 

-0.273-4 


* 

-0.2072 

0.0025 

-0.1S32 

-0.8705 

0.0343 

9. 1357 

1.6220 

0.0455 



-0.6S27 

-0.1162 

0.0234 

-0.1819 

-0.5302 

-0.2734 

0.0433 

1.CC63 


\ 
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Finally, to investigate tho effects of process and measurement noi3e on 
the algorithm performance, two cases were run. The first case included 
white Gaussian noise (WGK) with a sigma of 0.1 (ft/sec* ) on the outputs 
only. The results are shown in Flguro 4-24. Littlo performance degradation 
was noted, so the sigma was increased to 0.25 (ft/sec 1 ) and process noise 
with sigma 0.1° was added to tho control feedback commands. The results of 
this tost are shown in Figures 4-25 and 4-26. The 'measured' outputs on the 
right in Figure 4-25 contain the measurement noise, while the outputs on the 
left are the actual vibrations prior to the addition of the measurement 
noise. Similarly, in Figure 4-26, the control tine histories on the left 
are the feedback commands, and those on the right are the controls actually 
applied to tho simulation, i.e. with tho process noi3e added. 
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Figure M- 13 . Power Spectra of the Steady-State Uncontrolled and Adaptively 
Controlled Accelerometer Outputs at 120 kts/1000 ft. 
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Figure k-18. Selected Trajectory Parameter Time Histories for Adaptive 
Control at 80 kts/1000 ft with Time-varying Rotor Angular 
Velocity 

























figure M-2Q. Selec'.ed Parameter Estimate Time Histories for Adaptive 

Control at 80 kt3/’000 ft with Tlme~varying Rotor Angular 
Velocity 
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Figure H-25. Control Inputs for AdaptLve Control at 120 kts/1000 ft 
with Measurement Noi3e (o ~ 0.25 ft/sec 2 ) and Process 
Hoise (o » 0.1°) Included 
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Figure 1 )“26. Accelerometer Outputs for Adaptive Control at 120 kts/1000 ft 
with Measurement Noise (cr - 0.25 ft/sec 2 ) and Process Noise 
(o » 0.1°) Included 
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4.3*3 Discussion of Adaptive Controller Simulation Results 

The simulation results of the adaptive controller applied to the GENHEL 
model of the RSRA helicopter indicate significant vibration reduction at all 
steady-state flight conditions tested. Significant vibration reduction was 
also achieved during changing flight conditions. Notably, the reduction was 
achieved with no unique a priori information about the helicopter or flight 
conditions. The algorithm Initialization and tuning were Independent of the 
flight conditions . 

Adaptive controller performance under steady-state flight conditions 
without measurement noise or process noise Included is shown in Figures 4—5 
through 4—17 . As discussed previously, the first second of the acceleration 
time histories indicates the uncontrolled, unforced vibration levels; from 
one to two seconds wide band noise is being added to the controls, but no 
feedback for vibration control is present; from two seconds on the control 
feedback is present. As can be seen by noting the times at which the vibra- 
tion levels drop suddenly to their final values, the time from control 
feedback initiation to parameter identification disabling (and consequently 
wide band noise turn-cff) is on the order of two to three seconds. This is, 
of course, a strong function of the ID-OFF threshold level. In Figure U-16, 
the results for the flight condition 120 knots at 5000 feet • indicate a tran- 
sient time in excess of 4 seconds. For this flight condition, the ID-OFF 
threshold was lowered from 3.0 ft/sec 2 to 2.0 ft/sec 2 . 

For all the steady-state flight conditions, less collective is used 
than swashplate control. It was also noted during the algorithm tuning 
phase of the analysis, near the stability boundary large amounts of .collec- 
tive were used. Correspondingly larger amounts of swashplate deflection 
were used as well. Apparently, attempts to use more collective deflection 
to decrease the vertical vibration result in increasing amounts of power at 
N/rev cross-coupling into the lateral channels; thus requiring larger 
swashplate deflections. Incorporating cross-coupling terms in the MISO 
parameter identification algorithm might alleviate these problems. As far 
a3 the GENHEL simulation is concerned, however, it is seemingly necessary to 
keep the amount of main rotor collective used for vibration control small. 
Increasing collective deflections are indicative of impending instability. 


- 99 “ 



During each of the 3 teady- 3 tate flight condition test 3 , it was also 
noted that the RSRA simulation could not maintain the supposed trim condi- 
tion. In each case, non-zero body angular rates developed which after 
several seconds caused slow divergence of the velocity and altitude pro- 
files. The most common divergence was due to a small positive pitch rate 
and resulted in (initially) increasing altitude and decreasing velocity 
profiles. For the 120 knot/7000 feet flight condition, the divergence wa 3 
rapid however, and over the eight second test interval, the flight condition 
varied significantly. After initially suppressing the vibration below the 
ID-OFF threshold, the vibration level increased to the point where the lo oh 
threshold was surpassed. Increasing amounts of control authority were re- 
quired throughout the eight second test. To prevent this divergence, the 
roll and pitch rates of the body were set to zero in order to prevent the 
divergence. Note this was required for the 120 knot/7000 feet flight condi- 
tion only. 

Table H-2 summarizes the final steady-state vibration and control 
amplitudes for the various flight conditions without measurement or prccee3 
noise. For the 3 teady state conditions, the Input amplitudes gradually in- 
crease for more demanding flight regimes, though not nearly as fast sas for 
the f>xed-gain controllers (cf. Table 3 — 1 ) . These results are highly de- 
pendent upon the ID-0N/0FF thresholds and the time history of vlae band 
noise input for robust parameter identification. This dependence is pos'- 
sibly a consequence of overpararaeterization/misinodeling as indicated by the 
large parameter correlations in Table h— 1 . 

The re 3 ult 3 of the adaptive control with time-varying n are shown in 
Figures 4-18 through U-20. The rather abrupt change in R results In a sig- 
nificant decrease in velocity over the interval ( 6 c was increased in an 
attempt to maintain altitude) as well. As can been seen in the output time 
histories, during the period of R-variation, the vibration suppression is 
less than impressive. This is a consequence of the time constants In the 
narrowband filters In the identifier and controller structures. The center 
frequency of the filters tracks the R-variatior. 3 , but the bandwidth regains 
fixed. The time constants are too large for the filters to be effective 
when the center frequencies are changing rapidlv After the fl-variatich 
ceases, the filters reach steady-state and the controller performance im- 
proves substantially. The bandwidth of the filters could certainly be 
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TABLE 4-2. SUMMARY OF ADAPTIVE CONTROLLER STEADY-STATE ACCELERATION AND 
CONTROL AMPLITUDES FOR STEADY-STATE FLIGHT CONDITIONS 


f Flight Condition 

Steady-State 

Steady-State Feedback Control 

VEQ 

Altitude 

RSS ACC 


Peak-to-Peak ] 

(knot3) 

(feet) 

(ft/sec 1 ) 

B l3 (deg) 

9 c (deg) 

A l3 (deg) 








100 

0.8 

0.2 

0.3 

0.7 

5 







1000 

1.0 

0.5 

0.3 

0.7 


100 

1.0 

0.7 

0.6 

1.8 

40 







1000 

1.0 

1.2 

0.4 

1 .6 


100 

1.2 

0.3 

0.4 

1 .5 


1000 

1 .0 

1 .0 

0.1 

0.4 

120 

2000 

1.0 

1.0 

0.2 

0.6 


5000 

1 .2 

0.9 

0.3 

1.1 


7000 

1 .5 

1.0 

.0.6 

1.6 


increased somewhat during period of significant R-variation. However, more 
sophisticated logic would be required, and caution would have to be exer- 
cised so as net to excite the nonlinearities of the helicopter, thus 
Invalidating the approximate model of the vibration dynamics. 

The vibration results during transition frem 80 to 125 knots are inter- 
esting (cf. Figures 4-21 through 4-23)). After the identification start-up 
procedure, the vibration is suppressed by a factor of 5“10 at 6 seconds when 
the maneuver begins. The controller is continuously operating and iden- 
tification comes on when the vibration performance measure exceeds the ID-CN 
threshold about one tenth cT a second after the start of the maneuver. 
There is a fair amount of dynamics due to pilot motion in the period from 6 
to 1C seconds. By viewing parameter plots in Figure 4-22, it is evident 
that at during the maneuver, the model poles become more heavily damped (a^ 
and a 2 decrease) and that the phase of the second-order transfer function is 
changing radically (b^ and b 1 change symmetrically by significant amounts). 
The vibration level during this time is slightly less than the uncontrolled 
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vibration, even with the additive excitation. Note, the pilot inputs during 
this period are severe enough that they have some contribution in the 
bandwidth around N/rev and thereby "fight" the vibration suppression con- 
troller. With more stable rigid body motion after ten seconds the 
parameters begin to stabilize and by sixteen seconds the vibration threshold 
is passed such that the identification turn3 off. As the control gains reach 
their steady-state values, the vibration level continues to decrease. Other 
small increases are due to some lateral acceleration at twenty and thirty- 
one seconds. The ID-ON threshold is not surpassed for the remainder of the 
test. 

The results of the adaptive controller tests with measurement noise 
only (Figure U— 24 ) and with measurement and process noise (Figures 4—25 and 
4 - 26 ) indicate that the algorithm performance does not degrade significantly 
with moderate amounts of noise added. Certainly, with broadband measurement 
noise only, little degradation is to be expected since the narrowband con- 
troller filters output most of the noi3o before generating tho feedback 
signals. Process noise clearly i3 the dominant factor (unless powerful nar- 
rowband measurement noise is present), and for wideband noise with c-0.1°, 
the performance degradation is minimal. 



SECTION 5 

SUMMARY OF FIXED-GAIN AND ADAPTIVE CONTROLLER RESULTS 


The design of a fixed-gain controller for suppressing vibration in the 
model (GENHEL) of the RSRA helicopter was successfully accomplished using 
frequency 3haped cost functionals and a linearized model of the vibration 
dynamics obtained by numerical calculation of stability derivatives in a 
fourteenth-order model of the vibration dynamics. The frequency shaping was 
necessary in order to prevent feedback at frequencies other than R (the 
region of model validity) which could cause instabilities. This controller 
was found to be destabilizing, however, when flexible fuselage modes were 
added to the simulation, resulting in the need for the inclusion of addi- 
tional modes in the model. Noting that the model had grown too complex for 
control of nearly periodic disturbances, an approximate second-order 
autoregressive with exogenous input (ARX) model of the input-output 
relationships was identified instead. Based upon thi3 approximate model of 
the vibration dynamics, a fixed-gain controller was designed at a flight co- 
ndition of 120 knots forward speed and 100 feet altitude. The performance 
of the controller at and away from the design point was quite good. 

The design of the adaptive controller employed a similar / RX model of 
the vibration dynamics and used the frequency shaping concept in the control 
design as well. The Implementation was signlf lcantly-mcre complex than the 
fixed-gain- controller, involving on-line parameter identification, state es- 
timation, and control gain calculation. The algorithm performed quite well 
at all flight conditions as well as during transients. Robustness to 
measurement and process noise was demonstrated as well. 

Though the temptation to compare the results of the two algorithms is 
great, this should be done in the light of several important polr-t3. First, 
and most important is the fact that the underlying assumptions upon which 
the two algorithms are based are differen t. The fixed-gain controller is 
designed assuming that the model of the vibration dynamics is fixed and 
known for all future times, whereas the adaptive algorithm a3saT.es the model 
form is fixed, but that the parameters in the model may vary with time. 
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This subtle difference is similar to the difference between causal and anti- 
causal filters. Knowledge of the future can always be used to great 
advantage in control design, and the fixed-gain controller is using such 
knowledge. The adaptive algorithm, on the other hand, attempts to identify 
the present system, estimate the current states, and control the future out- 
puts based on the past and limited expectations concerning the future. 

A second key difference in the controllers involved the 3et of con- 
trols used for vibration suppression. All four controls were used in the 
fixed-gain design- whereas . only three were used in the adaptive algorithm; 
tall collective wa3 not used. Another important difference wa3 that dif- 
ferent models were used for system identification. The fixed-gain control 
design was based on a fully HIM3 identified ARX model structure. The adap- 
tive controller ignored the output coupling identifying parameters in three 
decoupled MISO ARX models. 

However, in spite of all these important differences, there are a few 
general comments which can be made. First, the algorithms both performed 
equally well at all the flight conditions tested in terms of overall vibra- 
tion suppression. Away from the design point of 120 knots/100 feet, the 
fixed-gain controller required more control authority than did the adaptive 
algorithm which is to be expected. However, the fixed-gain controller ex- 
hibited faster convergence than the adaptive controller at and near the 120 
knots/luO feet design point. This is a direct consequence of the different 
assumptions discussed above.. 

As a final note, it should be remembered that both controllers were 
controlling a simulation and not a real helicopter. The nonlinearities of a 
real helicopter are most assuredly more severe than those in the GENHEL 
simulation of the RSRA, and real actuator dynamics (not modeled in the 
simulation) will certainly have an effect on the controller performance. 
However, the most important point with respect to unmodeled nonlinearities 
is that the objective is to control vibration at or near a few select fre- 
quencies . only. Global validity of the models employed is not required. 
Approximate models, valid in the frequency ranges of interest, are suffi- 
cient for effective control designs, both fixed-gain and adaptive. 
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SECTION 6 
RECOMMENDATIONS 

The results of the adaptive controller tests are quite promising, 
however further analysis Is warranted before stronger conclusions can be 
reached. A more extensive process and measurement noise sensitivity 
analysis should be conducted in order to ascertain the performance degrada- 
tion with increasing measurement noise. Since the adaptive algorithm 13 
nonlinear, there is little that can be done analytically; Monte Carlo 
analysis is appropriate for these investigations. 

Further trade-off studies should be performed in which the effect of 
the bandwidth of the noise added for identification on the transient vibra- 
tion levels is ascertained. Noise outside the controller bandwidth (which 
is significantly smaller than the noise filter bandwidth) is not controlled 
and adds to the vibration level during the ID process. On the other hand, 
as the noi3e filter bandwidth decreases, the identif lability of the model, 
parameters decreases resulting in larger variance estimates and possible 
subsequent controller degradation. 

As discussed in Section 4, the controller implemented used the fil- 
tered state estimates in the feedback loop. For real-time implementation, 
this will resul in a delay in the control actuation due to the finite com- 
putation time. Use of the predicted state estimates eliminates the problem 
of the delay, but introduces states with a larger variance in the control 
loop. The effects of the delay versus the predicted state feedback should 
be Investigated. The fundamental factor in determining which approach is 
superior is the adequacy of the second-order model of the helicopter 
dynamics around N/rev. If the model is not good, the predicted states will 
deviate significantly from the 'true' states and the controller performance 
will degrade. 

Finally, further studies should be performed on a simulation including 
nigher order aeroelastic nonlinear effects and, more importantly, with more 
realistic models of the actuators. It is anticipated that the later will 
have a more significant impact on the achievable performance than the 
former. Regardless of the magnitude of the nonlinearities, the second-order 
model employed is sufficiently accurate in the frequency ranges of interest 
for periodic disturbance suppression. However, control actuator lags and 
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limits could potentially limit severely the performance of both the fixed- 
gain and adaptive controllers. 

In spite of the caveats discussed above, the autotuning type of adap- 
tive controller successfully demonstrated in this study Is quite attractive 
for an operational environment, where manufacturers will not want to do 
identification on each vehicle that is produced or 'retune' the controller 
during periodic maintenance to track changes in the flexibility effects of 
structual members. Thus, both the fixed-gain and adaptive frequency shaped 
control techniques analyzed here should be tested in wind tunnels. A logi- 
cal sequence would be: 

Wind tunnel testing 

1) with a fixed-gain scheduled control law, 

2) with a single-point robust fixed-gain control law, and 

3) with an adaptive autotuning control law. 

Should the te3ts prove successful, demonstrating sufficient robustness to 
disturbances of all types, flight tests would be warranted. A similar se- 
quence a3 that outlined above would be appropriate. 


REFERENCES 


[1] Kalctka, J. t .'Rotorcraft Identification Experience', AGARD Lecture 
Series No. 104, London, November 1979. 

[23 McCloud, J.L. , III, and Kretz, M., 'Multicyclic Jet-Flap Control for 
Alleviation of Helicopter Blade Stresses and Fuselage Vibration' , 
Rotorcraft Dynamics, NASA SP-352, 1974, pp. 233-238. 

[33 Johnson, W. , 'Self-Tuning Regulators for Multi-Cyclic Control of Heli- 
copter Vibration', NASA Technical Paper 1996, March 1982. 

[4] Johnson, W., Helicopter Theory . Princeton University Press, 
Princeton, New Jersey, 1980. 

[5] Gupta, N.K., 'Comparison of Frequency-Domain and Time-Domain Rotorcraft 
Vibration Control Methods' , NASA Contractor Report under Contract NAS2- 
11271 , Apr. 1984. 

[6] Gupta, N.K., and DuVal, R.W., 'A New Approach to Active Control of 
Rotorcraft Vibrations', AIAA Journal of Guidance and Control , Nov. -Dec. 
1981. 

[7] Gupta, N.K., 'Frequency-Shaping of Cost Functionals: Extensions of LOG 

Design Methods', AIAA Journal of Guidance and Control , Vol. 3, No. 6, 
Nov. -Dec. 1980, pp. 529-535. 

[8] Gregory, C.Z., Jr., and Gupta, N.K., 'Rotorcraft Vibration Controller 
Design and Evaluation,' ISI Report No. 14, ’repared by Integrated 
Systems, Inc. for NASA-Ames Research Center under Contract NAS2-10951, 
April 1982. 

[93 Bierman, G.J., Factorization Methods for Discrete Sequential 
Estimation , Academic Press, New York 1977 

[103 Bryson, A.E., and Ho, Y.C., Applied Optimal Control , Hemisphere 
Publishing Co., Hew York, 1975 

[113 Franklin, G.F., and Powell, J.D., Digital Control of Dynamic Systems , 
Addison Wesley, Menlo Park, California, 1980 

[123 Doyle, J. C. and Stein, G., 'Multivariable Feedback Design: Concepts 
for a Classical/Modern Synthesis', IEEE Trans. Automat. Contr., AC-26, 
pp 4-16, February 1981. 

[133 Safanov, M. G., Laub, A. J. and Hartmann, G. L., 'Feedback Properties 
of Multivariable Systems: Ihe Role and U3e of the Return Difference 

Matrix', IEEE Trans. Automat. Contr., AC-26, pp 47-65, February 1931. 

[143 Kosut; R. L., 'Robustness Issues in Adaptive Control', IFAC Workship on 
Adaptive Systems in Control and Signal Processing, San Francisco, CA, 
June 20-22, 1983. 


107 - 



[15] Molusia, J .A . .Hammond , C.E., and Cline, J.H.,’A Unified Approach to the 
Optimal Design of Adaptive and Gain-Scheduled Controllers to Achieve 
Minimum Helicopter Rotor Vibration', 37th Annual Forum of the American 
Helicopter Society, New Orleans, La., May 1981. 



APPENDIX A 

N ASTRA M MODAL ANALYSIS OF RSP.A FLEXIBLE FUSELAGE 




qUaUi 

<n£ pcu.. 


The first 34 natural frequencies, generalized masses and generalized 


stiffnesses for the 

RSRA flexible fuselage 

are 

shewn In Table A-1. -The 

eigenvector of mode number 15 , tl\ 

j closest mode 

to 

the rotor frequency H/rev 

is shown in Table A-2 

. In Table 

A-2, point ID 

number 1 1 Is the 

rotor frame 

(HUB) 

and point number 546 is considered to be the 

pilot location 

• 


Table A-1 . 

MODAL ANALYSIS OF THE RSRA 

HELICOPTER MODEL 




(MARCH 4, 

1981 - NASTRAN 

12/16/81 ) 



REAL 

EIGENVALUES 




»!w>r 

EIGENVALUE 

RADIANS 

CYCLES 


GENERALIZED 

GENERALIZED 

'0. . 





HAj; 

stiffness 

1 

0.0 

0.0 

0.0 


R.:5*' e 2E40l 

0.0 

•) 

0.0 

0.0 

0.0 


4.o54'92E40‘. 

0.0 

3 

0.0 

0.0 

0.0 


4.523220E401 

0.0 • 

4 

0.0 

0.0 

0.0 


1.778937E401 

0.0 

5 

0.0 

0.0 

0.0 


1.921 130E402 

0.0 

6 

0.0 

0.0 

0.0 


1.293919E402 

0.0 

7 

7.621417E402 

2.74072ES401 

4.393035E400 


1.42221CE400 

234393403 

B 

9.63015GE402 

3.103250E401 

4.9239/5E400 


1.267294E400 

1.220424E403 

9 

2.3083981403 

4.B04579E401 

7.444724E400 


1.712277E400 

3.952422E403 

10 

3.559379E403 

5.944053E401 

9.49524CE4CO 


1.545444E+00 

5.572724E403 

11 

4.190244E403 

4.473225E401 

1.030244E401 


2.872793E400 

S.203774E404 

12 

6.S77077E+03 

8.10991EE401 

1.290734E401 


1.B5S145E400 

1.222129E404 

11 

8.733301E403 

9.348154E401 

1.487805E401 


1.67215SE400 

1.434043404 

14 

9.494004E403 

9.743719E401 

1.550741Et01 


2.87342EE400 

2.728034E404 

15 

1.230904E404 

1.109441E402 

1.745742E401 


3.987P81E400 

4.902542E404 

14 

1.369259E404 

1.170154E402 

1, 84235 TEfOi 


2.864976Et00 

3.925434E404 

17 

1.448761E404 

1.203445E402 

1.915440E401 


3.575121E4C0 

5.179495E404 

IS 

1.591GS2E4C4 

1.241579E402 

2.007S44E401 


4.220672E+00 

4.717550F404 

19 

1.732225E404 

1.3I4140E402 

2.094702E401 


3.140145E400 

5.474117E404 

20 

1.9B01C4E404 

1.407142E402 

2.239543401 


4.230754E400 

8.377333E404 

21 

2.452744E404 

1.544125E402 

2.492545E401 


2-129577E+00 

5.22331 1E404 

m 

2.753163E404 

1.459244E402 

2.440904E401 


3.311485E400 

9.117037E404 

23 

2.805484E404 

1 . 675O1EE+02 

2.445375E401 


1.987014E400 

5.574944E404 

24 

2.E83033E404 

1.499421E402 

2.704713E401 


2.27747SE400 

4.57S009E404 

25 

3.316444E404 

1.821114E402 

2.893394E401 


3.017274E4CO 

1.000449E4C5 

26 

3.497404E404 

1.870U5E4Q2 

2.974412E401 


3.317150E40.' 

1.160141E40G 

27 

3.473318E404 

1.914590E402 

3.050348E401 


5.22SC42EHV' 

1.920424E405 

28 

*.007387E404 

2.001846E402 

3.1P4037E401 


2.8187ME40'- 

1.1295B7E405 

29 

4.709474EM4 

2.1701789402 

3.453946E40! 


3.761E31E40D 

1.771700E405 

30 

5.05396SE404 

2.2481035402 

3.57794SE401 


3.510554E407 

1.7742I4E405 

31 

3.095703E404 

2.257344E402 

3.592710E1C1 


5.241S22E400 

2.471077E405 

32 

5.438099E404 

2.374468E402 

3.779033E401 


3.847402E400 

2.149204E405 

33 

5.498542E404 

2.3S7142E+02 

3.799234E401 


4.73177CE400 

2.496419E+05 

34 

5.939238E404 

2.437045E402 

3.878710E401 


8.498703E400 

5.164411E405 
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Original p'ag^ is: 
. QE poor QUALITY 


Table A-2. EIGENVECTOR OF MODE 15 OF RSRA HELICOPTER MODEL 
(MARCH 11, 1981 - NASTRAH 12/16/81) 


EIGENVALUE* 1.23J904EH4 

CYCLES * 1.765’42E401 REAL EIGENVECTOR HO. 15 

POIWT ID. Tl T2 T3 R1 

1 -2.18?307E-C2 -2. 3321 10E-02 -5.G30057E-03 -I.112256E-03 -2.3®9»E4E-G4 -4.e2i2'TE-05 

2 -2. 1041 B3E-C-2 -2.520769E-02 -1.45G342E-P2 -1.1I2256E-03 -2.30908SE-54 -4.4:42*7E-:-5 

3 -2.274430E-02 -2.520769E-02 2.44224CE-02 -1.U2256E-03 -2.399984E-04 -4.42627?E-05 

8 -3.012502E-02 1.3P9802E-02 -9.998965E-04 -1.4I5554E-03 -5.333254E-04 -5.634?°5£-05 

9 -2.431705E-02 -1.301S5SE-02 -9.9S3717E-04 -1.112256E-03 -2 .399724E-04 -4.424272E-05 

10 -2.452504E-02 -2.705824E-O3 -9.9S835IE-04 -1.I12256E-03 -2.399986E-04 -4.624277E-05 

11 -4.678098E-02 3.128137E-02 -1.001022E-03 -2.032791F-03 -1.13C329E-03 -4.554213E-05 

12 -2.454072E-02 6.748619E-03 4.1144D1E-02 -9.S51572E-04 B.1I2236E-04 -1.45S077E-04 

13 -2.812449E-02 7.145543E-03 4.O94339E-02 -1.475374E-03 A.445037E-04 -1.405402E-04 

101 -4.249715E-03 1.41C412E-01 -1.449113E-01 -2.MG294F-02 -2.720051E-03 -I.054095E-03 

527 -2.142647E-02 5.115342E-02 4.02201CE-02 -2.393977E-03 -4.4435HE-04' 0.0' 

535 2. 2279 ICE-02 4.994913E-02 -7.478<43E-02 -1.340232E-03 -1.6-57347E-03 0.0 

533 2.210339E-02 3.473072E-02 -3.324542E-02 -1.43I0I7E-O3 -1.53625CE-03 0.0 

546 -2.U3531E-02 3.709319E-02 5.521093E-02 -1.D37675E-03 -5. 1T.G593E-04 0.0 

447 2.1950672-02 2.6337002-02 -4.299799E-C2 -6.644740E-04 -1.331C93E-03 0.0 

650 -2.0S3491E-02 2.998214E-02 5.895294E-02 -1.543360E-03 -5.3121292-04 0.0 

KOI' 1.147524E-02 -2.G48305E-02 -3.601325E-C3 0.0 -3.534379E-04 0.0 

1402 1.112423E-02 -3.05451C2-02 -4.8333ME-03 0.0 2.CS8741E-04 0.0 



1403 

1.150517E-02 

-3.212431E-02 

-4.830200E-03 

0.0 

2.776334E-04 

0.0 


. 1404 

1.25S090E-02 

-3.333530E-02 

-4.G21813E-03 

0.0 

-1.319432E-04 

0.0 


1440 

1.761409E-02 

-2.337344E-02 

1.162864E-02 

0.0 

3.470051E-04 

0.0 


1441 

1.4U692E-02 

-2.462533E-02 

5.954643E-03 

c.o 

4.553997E-04 

0.0 


1442 

1.229327E-02 

-2.626352E-02 

5.02792'E-M 

0.0 

3.68394BE-04 

0.0 


1504 

1.514671E-02 

-3.134985E-02 

-5. 164 4592-03 

0.0 

3.024370E-O4 

0.0 


154G 

2.324402S-O2 

-2.772009E-02 

8.201564E-03 

C.O 

-3.217202E-04 

o.o 

! 

1601 

1.979644E-02 

-2.5251382-02 

-8.438348E-04 

0.0 

J.935107E-O5 

0.0 


1602 

1.846534E-02 

-2.53077GE-02 

-5.22796DE-03 

0.0 

4.962934E-04 

0.0 


1603 

1.760415E-02 

-2.551403E-02 

-9.697iB5E-03 

0.0 

4.B54839E-04 

0.0 


1604 

1.631994E-C2 

-2.60851EE-02 

-1.640774E-02 

0.0 • 

6.136770E-04 

0.0 

- 

1640 

2.307229E-02 

-2.9294S0E-02 

2.410931E-02 

0.0 

-8.V24071E-04 

0.0 


1641 

2.142116E-02 

-2.718537E-02 

1.093953E-02 

0.0 

5.3S8230E-04 

0.0 


1642 

2.043333E-02 

-2.603326E-02 

4.281957E-03 

0.0 

5.32322BE-G4 

0.0 


4913 

6.149341E-02 

8.3176A4E-01 

1.570026E-0I 

0.0 

0.0 

0.0 

i * 

5101 

-4.3E4285Et00 

1.724457E+00 

-5.274023EtC0 

0.991367E-02 

M43254E-03 

-6.650239E-02 

» 

5102 

-4.586285E400 

1.164374E400 

-5.303935E+00 

8.991367E-02 

4.145254E-03 

-6.65C239E-02 


5103 

-4.5S6285E400 

2.1CB424E-01 

-5.3483A9E+00 

E.991367E-02 

4.H5254E-03 

-6.650239E-02 

\ 

5130 

4.875£5iE100 

2.01 4762E-0I 

5.531722E+00 

9.158336E-C2 

-i.52754I.E-C3 

-6.7C.'73;E-02 

j 

5131 

4 . S75351 E+OG 

.1.1625152+00 

5.309815E+00 

9.155336E-02 

-1 .527947E-G3 

-O.7C2728E-02 

f 

i 

5132 

4.875351E100 

1.727019E+00 

5.494946E+00 

9.133336E-02 

-1.527945E-03 

6.70273BE-02 
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APPENDIX B 

LINEARIZED MODEL OF THE RSRA VIBRATION DYNAMICS 






Thi3 appendix presents the dynamic matrices for the identified linear 
model of the RSRA vibration dynamics in steady-state flight at a l2u-knot 
forward speed and 100 feet altitude flight condition. The discrete domain 
MIMO ARX model given in equation 3.1; 


'k+i 


- Vk * 


Vk-1 


B u 
ik 


+ B. u, , 
i k~1 


♦ r w, 


(3.1) 


where u is the vector of control inputs, y is the vector of outputs and 
w is a vector of the sine and cosine of the rotor angle (i l>) , was used to 
perform least-squares identification of the coefficient matrices. 
Transforming this equation into state-space form for controller design was 
performed. In detail in Section 3 (cf. equations 3.21 through 3.23). Using 
the notation employed in Section 3.2, the following matrices represent the 
identified model of the RSRA vibration dynamics (i.e. the plant) coupled 
with the freqr.ency shaping filter dynamics. 


i 

; 

• i 




,:.-J 


F - 


L-: 


CO. IKS 1 T1RU 6 


1.53405403 

0.4314D-03 

2.34105-03 

1.0500D400 

0.00005400 

0.05005400 

-7.73075-03 

1.93155400 

2. 10245-03 

0.00005400 

1.00005400 

0.0000P400 

2.47535-04 

-5.23005-02 

1.97905400 

0.00005400 

0.00005400 

1.00005400 

-9.91 395-01 

-9.97475-03 

-2.44935-03 

0.00005400 

0.00003400 

0.00005400 

0.45325-03 

-9.9530&-01 

-2.57315-03 

0.00005400 

0.00005400 

0.00005400 

2.87545-03 

4.31415-02 

-9.90495-01 

0.00003400 

0.00003405 

0.00003400 

1.21445-02 

5.23425-05 

1.43325-05 

6.12215-03 

0.00005400 

0.00005400 

2.42475401 

1.05385-01 

2. 04355-02 

1.22325401 

0.00005400 

0.00005400 

H. 73235-05 

1. 21315-02 

1.23715-05 

0.00005400 

4.12215-03 

0.00005400 

-9.45405-02 

2.42335401 

2.57145-02 

0.00005400 

1.22325401 

0.00900400 

1.51575-04 

-3.23245-04 

1. 21165-02 

0.05005400 

0.00005400 

6.12215-03 

. 3.02335-03 

-4.45335-01 

2.42075401 

0.00005400 

0.0.0005400 

1.22323401 

CCU5WS 7 TIOTiU 12 
0.00005400 0.00005400 

0.00005400 

0.CWD400 

0.00005400 

0.50002400' 

0.05505409 

0.00005400 

0.00055400 

0.00005400 

0.05000400 

0.00095400 

0.03005400 

0.C0I054C0 

O.COO0D40O 

0.00000400 

0.00003400 

0.00005400 

0.00005400 

0,00005403 

0.00005400 

0.00003400 

0.030SD400 

0.00505400 

0.00055400 

0. 00005400 

0.00003400 

0.00500400 

0.00035 490 

0.00005400 

0.000*5400 

0.00035400 

0.C053D4C0 

0.05555100 

0,05003400 

0.00003400 

9.93533-01 

9.97945-04 

0.00055403 

0.C-05334C0 

0.00003400 

0.00005400 

-1.22323401 

9.93885-01 

0. 000034 00 

0.005S3400 

0.0005D400 

0.00503400 

0.00005400 

0.00005400 

9.93325-01 

9.97763-04 

0.00003400 

0.00055400 

0.00005400 

0.00005400 

■ 1.22325401 

9.93SC3-01 

0.00005400 

0.03595400 

0.00003400 

0.00095400 

0.05005400 

0.05005400 

9.92885-01 

9. 97962-04 

0.00005400 

0.00595450 

0.00005400 

0.00093400 

-1.22223401 

9.92235-01. 


-1.37375-02 
7.19735-02 
2.00295-01 
9.00fcIH)3 
-7.24105-02 
-2.14415-01 
-3.41955-05 
-1.63221. -01 
4. <0435-04 
0.C9J5D-01 
1.22625-03 
2.44995400 


-2.22175-01 
-1. 2555-01 
-1. 91275-01 
2.13225-01 
1.08755-01 
1.53715-01 
-1.3443D-03 
-2,72935400 
-4.44575-04 
-1.32735:00 
-1.17115-03 
-2.33935400 


1.90155-01 
4.49025-02 
9.43145-02 
-1.B754D-01 
-3. 90555-02 
-B. 02235-02 
1. 14415-03 
2.32593400 
2.B714D-04 
5.73495-01 
5.30475-04 
1.15935400 


-4.13455-03 
-1.21165-02 
2.50115-02 
4.07335-03 
1. 20155-02 
..91145-02 
-2.55305-05 
-5. 12035-02 
-7.41855-05 
-1.43225-01 
1.53125-04 
3. 05935-01 


I 

I / I 
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This model wa 3 U3ed for control de3lgn as discussed in Section 3 . 2 . The 
control gains used and the open-loop and closed-loop eigenvalues of the 
linearized vibration dynamics ares 


Open-Loop Eigenvalues 


0 

U 

w 

n 

C 

( 1 /sec) 

(rad/s) 

(Hz) 


4.25440*00 

1.24610+02 

1 .98340+01 

3.41 40D-02 

2.57470+00 

0.7334D*O1 

1.39130*01 

2.77200-02 

4.40220+00 

1.00110+02 

1.59560*01 

*.37325-02 

4.25460*00 

-1 .24610+02 

i.m<D+oi 

3. *1*00-02 

2.57670+OO 

-8.733*0+01 

1. 37130*01 

2.97201'-©: 

4.4022D+00 

-1.00110+02 

1.37540*01 

*. 3732I'-02 

-1.02380-0? 

5.10710*02 

1.74270*01 

-7.24730-12 

5. 07030-12 

1.1071 no: 

1. 76270*01 

5.302*0-14 

3.09270-11 

1 . 10710*02 

1. 74270*01 

3.51410-13 

-1.02325-0? 

-1.10710+02 

1.76270+01 

-9.2417D-12 

5.76740-12 

-1. 10710*02 

1.74270*01 

5.3902D-14 

3.C302D-11 

-1.10712102 

1.74270+O1 

3.50*90-13 


Closed-Loop Eigenvalues 


0 

( 1 /sec) 

0 ) 

(rad/s) 

hi 

n 

(Hz) 

C 

7.1*153*00 

1.24405*02 

2.01595*01 

5.64110-02 

7.14150100 

-1.26400*02 

2.0157D+01 

5.6111P-02 

*.*5340*03 

0.65723*01 

1.33040101 

5 .13760—02 

*.*5340*00 

-3.45725*01 

1.380*0*01 

5.13760-02 

6.67555tW 

9.034201 01 

1.5496C101 

6.772*0-02 

6.675SB+M 

-9.03*20+01 

1.56740+01 

6.772*0-02 

5.902*0*00 

1.12400102 

1.7923D101 

5.2**00-02 

5.902*0*00 

-1.12403+02 

1.77230+01 

5.2**00-02 

*.97*30*00 

1.09470*02 

1.7*810*01 

*.5*730-02 

*.99*30*00 

-1.07470*02 

1.7*810*01 

*.5*730-02 

1.02500+00 

i. 10800*02 

1.76*40*01 

1.6*690-02 

1.B250D+00 

•1.10500+02 

1.74140+01 

1.6*690-02 
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I 

i 

» 

i 

i 


i 

> 

i 



cousss i Ttsu t 

5.23920-02 3.94590-31 -1.2B31D-01 
-2.55S0D-01 6.24256-01 -1.G679DH31 
-4.3321&-02 5.00040-01 -2.S159D-02 
5.0362D— 02 -t.4042P-02 -3.05476-02 

CCLUtKS 7 TR.1U 12 

-3.03766-04 -4.94250-04 1.52236-02 
3.E0360-03 5. 91480-04 -9.1947D-03 
-0.90000-03 4.0:306-05 -6.03026-03 
-1.42166-03 -1.91956-05 3.91016-03 


7.63336-02 3.22056-01 -1.42016-01 
-2.39416-61 7.07906-01 -1.90976-01 
-4.37410-02 4.V9430-01 -2.24370-02 
5.73126-02 2.02600-03 -3.59990-02 


7.35130-05 3.C940D-03 2.54156-05 ' 

-1.15790-04 2.22' 60-03 3.66100-05 

3.63396-06 1.13635-03 9.76650-06 

-1.95696-05 -2.25520-05 1.10470-05 . 




APPENDIX C 

DIRECT FEEDTHROUGH MODIFICATION TO LQR COST FUNCTIONAL 


In this appendix, a derivation is presented indicating the modifica- 
tions to the standard LQR c03t function in regulator design required to 
account for a direct feedthrough term, i.e. the presence of a component of 
the inputs in the outputs. The derivation is given in the continuous 
domain to simplify notation. Including the direct feedthrough term, the 
dynamical model of the system to be controlled can be written: 


x - Fx + Gu 
y - Hx + Du . 


(C..1) 


The co3t function to be minimized is the standard quadratic functional: 


2 | Cy T Ay + u T Buj dt 


Minimizing .J subject to the constraints of the dynamical model by the method 
of Lagrange multipliers lead3 to the deflntion of the Hamiltonian and the 
Euler-Lagrange equations: 


Defining: 


H - y T Ay + u T Bu + X T (Fx ♦ Gu) , 


3H 

3H 

IS 


X - - x T h‘aH - U T D T AH - X T F , 


T T t T T T 

0 - x H AD + u B + X G + u D AD 


(C.2) 

(C.3) 

(C.H) 


B 

A 

F 


B ♦ D*AI. , 
A 


adeT’d t a 


F - G3 _1 D T AH, 


equation C . h can be solved for u to give: 

— 1 T — 1 T 

u - - B G X - B D AHx . 


(C.5) 


(C.6) 
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Substitutes equation C .6 into C.4, the Euler-Lagrange equations can be 
written in the standard form: 


• “ 1 T 

X - Fx - GB Cl , 

X - - FX ~ H T AHx . 


(C.7 ) 


Thus, with the substitutions given in C.5, the LQR design with direct 
feedthrough of inputs to outputs can be cast in the standard cost functional 
form which which leads to the forward Ricatti equation and associated feed- 
back gains: 


S - SF * F T S - SGB - 1 0 T S «■ H T AH , 

C - - B - 1 G T S , (C. 8 ) 

u - Cx . 

In order to obtain the feedback appropriate for the original system in terms 
of the regulator gains computed in C. 8 , equation C .3 is rewritten as 
follows: 


u - - Cx - B ^D^AHx - - Cx , 

C - C + b -1 d t ah . 


(c.9) 


It should be emphasized that even In the case of time-invariant system 
matricies (F,G,H,D), the 3teady-3tate solution to the forward regulator 
problem will not be the same a3 the fixed-gain regulator solution discussed 
in Section 3 . They are solving two DIFFERENT problems! In the forward, or 
causal regulator solution, only knowledge of the past is a33umed. The 
fixed-gain regulator assumes knowledge of the system over the entire inter- 
val which is usually assumed to be infinite. 
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